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SECTION  I 
INTRODUCTION 

There  are  two  commonly  accepted  means  of  spreading  flame  through  a com- 
bustible mixture  in  practical  combustion  systems.  The  first  of  these,  laminar 
flame  propagation,  depends  on  heat  conduction  and  the  diffusion  of  chemically 
active  species  into  the  adjacent  fuel-air  mixture  to  propagate  the  fire.  One  foot 
per  second  is  a typical  laminar  flamespeed  for  stoichiometric  hydrocarbon-air 
mixtures.  The  second,  turbulent  flame  propagation,  adds  the  turbulent  transport 
of  small  elements  of  flame  a short  distance  into  the  unburned  mixture  to  act  as 
new  ignition  sources.  Turbulent  flamespeeds  in  hydrocarbon-air  mixtures 
typically  range  from  2 to  about  20  feel  per  second.  In  the  past,  when  attempts 
have  been  made  to  increase  the  turbulent  flamespeed  to  higher  values,  the  pres- 
sure drop  required  to  produce  the  turbulence  has  been  prohibitive  for  practical 
applications,  or  flame  stabilization  problems  prevented  operation  at  higher 
velocities.  Recently,  tests  in  a combustion  centrifuge  have  demonstrated  that 
centrifugal  force  can  be  used  to  increase  flame  propagation  rates  by  an  addi- 
tional factor  of  4 or  more.  This  report  describes  the  analytical  and  experi- 
mental program  conducted  under  Contract  F44620-73-C-0061  with  the  Air  Force 
Office  of  Scientific  Research  to  provide  a sound  basis  for  applying  this  phenom- 
enon to  Air  Force  combustion  problems. 
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SECTION  II 

RESEARCH  PROGRAM 

A 2-year  combined  analytical  and  experimental  effort  was  aimed  at  a thorough 
understanding  of  the  effects  of  centrifugal  force  on  combustion  and  the  development 
of  a mathematical  model  that  can  be  used  as  a design  tool  for  Air  Force  combustion 
systems.  The  first  year's  work  consisted  of  obtaining  experimental  data  with  the 
combustion  centrifuge  and  comparing  the  results  with  those  fi’om  other  sources.  A 
preliminary  mathematical  model  was  also  developed.  During  the  second  year 
further  experimental  data  on  the  behavior  of  buoyant  fluids  in  a dense  medium  were 
obtained  and  refinement  of  the  preliminary  model  was  completed. 

A.  EXPERIMENTAL  AND  ANALYTICAL  RESULTS 

The  experimental  and  analytical  work  consisted  of  four  phases:  (1)  the 
spreading  of  flame  through  a combustible  mixture  by  the  buoyant  movement  of 
flame  bubbles,  (2)  the  extinction  of  flames  in  high  centrifugal  force  fields,  (3)  the 
effects  of  Reynolds  No.  on  turbulent  flamespeed,  and  (4)  application  of  the  results. 
The  four  phases  will  be  discussed  separately. 

1.  Flamespreading 

When  a hydiocarbon  fuel-air  mixture  burns,  the  heat  released  by  the 
chemical  reaction  raises  the  temperature  of  the  burped  products  and  causes  them 
to  expand.  The  expansion  results  in  a lower  density  of  the  burned  products  than 
for  the  unburned  mixture.  Gravitational  and  centrifugal  forces  t^nd  to  organize 
fluids  of  different  densities,  with  the  heaviest  on  the  bottom  and  the  lightest  on 
top.  This  buoyant  force  causes  the  flame  to  rise  from  a burning  candle  in  'he 
earth's  gravitational  field.  In  more  intense  combustion  systems,  such  as  in  a 
turbine  engine  burner,  the  buoyant  forces  produced  by  the  earth's  gravitational 
field  are  insignificant  compared  to  the  forces  produced  by  the  flowing  air  stream 
and  do  not  measurably  affect  the  combustion  process.  If  the  earth's  gravita- 
tional field  is  supplemented  by  a centrifugal  force  field  several  thousand  times 
stronger,  the  buoyant  forces  again  become  predominant  and  completely  control 
the  flame  spreading  process.  The  resulting  buoyant  force  can  be  expressed  os: 

Fl>  * <'»>  - Tj)  Eq  (1) 
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where: 


Fg  is  buoyant  force  per  unit  volume 
a is  gravitational  or  centrifugal  acceleration 
is  density  of  cold  air 
is  density  of  hot  products  of  combustion. 

This  force  is  relatively  weak  in  the  earth's  gravitational  field.  For  example,  it 
amounts  to  approximately  0.  063  pound  for  a 1-foot  cube  of  hot  gas  at  ambient 
pressure  having  a temperature  of  3500°F  surrounded  by  air  at  80°F.  In  a 3000-g 
centrifugal  field,  however,  it  becomes  a force  of  190  pounds  on  a cubic  foot  of 
gas.  This  large  force  completely  dominates  the  movement  of  the  flame  within 
the  combustor. 

A series  of  basic  experiments  was  conducted  using  a combustion  centri- 
fuge. The  centrifuge  was  constructed  from  a 2. 65-inch  internal  diameter 
(nominal  3-inch,  heavy  wall)  stainless  steel  pipe,  6 feet  long  and  closed  at  both 
ends  with  weld  caps.  An  automotive-type  spark  plug  is  mounted  at  one  end  be- 
hind a turbulence  generator,  and  a 2000-psi  burst  disk  is  located  at  the  other 
end.  The  turbulence  generator,  consisting  of  a flat  disk  with  holes  in  it,  has 
40'/}  blockage  and  is  placed  2-1/2  inches  from  the  spark  plug  to  ensure  the 
generation  of  a fully  turbulent  flame  before  the  fire  reaches  the  instrumented 
section.  Ionization  probes  are  installed  at  the  locations  shown  in  figure  1,  and 
a pressure  transducer  and  a thermocouple  are  mounted  near  the  center  of  the 
pipe.  Flame  propagation  rate  measurements  are  made  between  stations  1 and 
11  to  (1)  provide  a long  lest  section  to  reduce  the  scatter  from  small,  local 
flame  speed  variations,  and  (2)  obtain  data  before  much  of  the  mass  is  burned 
and  the  pressure  rise  due  to  combustion  becomes  significant.  Slip  rings  permit 
the  continuous  recording  of  12  channels  of  information.  An  insulated  brass  disk 
and  roller  transmit  the  ignition  energy  from  a variable  energy  ignition  source  to 
the  spark  plug.  The  entire  unit  is  mounted  on  tapered  roller  bearings  and  can 
be  rotated  by  a 400-horsepower,  variable-speed  electric  motor.  The  assembly 
is  shown  in  figure  2. 
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on  Test  Stand 
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The  ionization  probes  are  steel  spikes  mounted  from  opposite  sides  of  the 
pipe,  with  one  spike  of  a pair  welded  to  the  pipe  and  the  other  inserted  through 
an  insulating  collar.  The  spikes  are  conical  in  shape  to  resist  bending  under 
the  high  centrifugal  loads,  and  the  points  are  separated  by  a 0.050-inch  gap  at 
the  center  of  the  pipe.  A 90-volt  electric  potential  is  applied  to  each  pair  of 
spikes,  and  the  flow  of  current  that  results  when  ionized  flame  completes  the 
circuit  is  recorded  on  an  oscillograph  at  a paper  speed  of  160  inches  per  second. 
Rotational  speed  is  recorded  on  the  oscillograph  and  also  indicated  on  a digital 
display  in  the  control  room.  The  signal  is  generated  by  redundant  reluctance 
pickups  mounted  adjacent  to  60-tooth  gears. 


The  fuel-air  mixture  is  prepared  in  a separate  tank  and  admitted  to  the 
centrifuge,  one  charge  at  a time.  The  mix  tank  and  centrifuge  are  first  evacuated 
to  below  50-micron  pressure  and  valved  off  to  check  for  leaks.  The  mix  tank  is 
then  isolated  from  the  centrifuge  and  fuel  is  admitted  to  a predetermined  pres- 
sure. Dry  compressed  air  is  added  until  the  desired  fuel-air  ratio  is  reached, 
and  an  internally  mounted  paddle  is  oscillated  to  stir  the  mixture.  One  full  lank 
can  supply  fuel-air  mixture  to  the  centrifuge  for  from  3 to  45  tests,  depending 
on  the  pressure  level  being  investigated.  The  centrifuge  is  enclosed  in  an 
insulated  housing  supplied  with  liquid  nitrogen  so  that  the  temperatures  of  both 
the  centrifuge  and  the  surrounding  gas  can  be  varied  from  +200oF  to  -100°F. 


To  conduct  a test,  the  centrifuge  is  evacuated  to  a 50-micron  pressure 
to  remove  the  products  of  combustion  from  prior  testing,  and  a fuel-air  charge 
is  admitted  to  the  desired  pressure  level.  The  mix  tank  is  then  valved  off,  and 
the  transfer  line  disconnected  and  capped.  The  pipe  is  then  rotated  to  the  de- 
sired speed  and  either  chilled  by  adding  liquid  nitrogen  to  the  enclosure  or  heated 
l>y  punning  it  in  the  enclosure  until  mechanical  energy  of  rotation  has  heated  both 
the  centrifuge  and  surrounding  air  to  the  desired  temperature.  Thermocouples 
in  the  enclosure  and  inside  the  centrifuge  monitor  the  temperatures.  After  the 
desired  temperature  is  reached,  the  spark  is  fired  and  the  data  are  recorded  on 
the  oscillograph.  The  results  for  propane-air  mixtures  are  presented  in  fig- 
ure 3.  The  observed  flame  propagation  rate  is  plotted  as  a function  of  the 
centrifugal  acceleration  at  the  igniter.  The  centrifugal  acceleration  midway  be- 
tween the  ionization  probes  is  about  two-thirds  of  this  value. 
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Figure  3.  Observed  Flamespeed  as  a Function 

of  Centrifugal  Acceleration  for  Propane 
Air  Mixtures 


FD  74050 


As  figure  3 shows,  thr  flame  propagation  rate  is  relatively  unaffected  by 
centrifugal  acceleration  up  to  about  200  g.  Above  this  value,  up  to  about  500  g, 
is  a transition  region.  At  values  above  500  g,  the  observed  flame  propagation 


rate  increases  about  in  proportion  to  the  square  root  of  the  centrifugal  force. 


A dashed  "square  root  curve"  is  added  to  figure  3 for  comparison.  Near  3500  g, 
a point  is  reached  where  the  observed  flamespeed  abruptly  reverses  the  trend. 
From  here,  it  decreases  with  increasing  centrifugal  acceleration  until  a limit  is 
reached  beyond  which  combustion  will  not  occur.  Fuel-air  mixture  ratios  other 
than  the  one  plotted  here  exhibited  the  same  behavior,  although  the  maximum 
flamespeed  and  flame  extinction  points  occurred  at  significantly  different  centrif- 
ugal acceleration  values. 


Tests  with  hydrogen-air  mixtures  yielded  somewhat  different  results.  A 
stoichiometric  hydrogen-air  mixture,  at  1-atm  precombustion  pressure,  showed 
a constant  flame  propagation  rate  over  the  range  of  centrifugal  acceleration 
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values  tested.  When  the  mixture  equivalence  ratio  was  dropped  to  0.6,  however, 
the  measured  flame  spreading  rates  corresponded  very  closely  to  those  for 
propane-air  mixtures  at  the  same  g values.  The  experimental  data  are  presented 
in  figure  4,  along  with  data  obtained  from  stoichiometric  hydrogen-air  mixtures 
at  0.33  *atm  precombustion  pressure.  Included  for  comparison  is  the  mean  curve 
through  the  propane-air  data  of  figure  3.  These  contrasting  sets  of  data  imply 
that  two  different  mechanisms  control  the  flame  spreading,  and  a hypothesis  is 
proposed  to  explain  them. 


Centrifugal  Acceleration  - g 


Figure  4.  Observed  Flamespeed  as  a Function  of  FD  74049 
Centrifugal  Acceleration  for  Hydrogen- 
Air  Mixtures 

Figure  5 shows  the  situations  schematically.  A flame  bubble  immersed  in 
a denser  fuel-air  mixture  (represented  by  the  cross-hatched  circles)  is  displaced 
upward  by  the  buoyant  force  resulting  from  a high  centrifugal-force  field  and  the 
density  difference.  In  an  increment  of  time,  At,  it  moves  a distance  equal  to  the 
product  of  the  time  increment  multiplied  by  the  velocity,  SB»  which  will  be  called 
the  "bubble  velocity."  During  the  same  lime,  turbulent  flame  propagation  has 
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caused  the  flamefront  to  advance  in  all  directions  by  an  amount  equal  to  the  pro- 
duct of  the  turbulent  flamespeed,  S^,  multiplied  by  the  same  time  increment. 
(Hot-gas  expansion  behind  the  flamefront  is  neglected  for  simplicity. ) The  ex- 
ample on  the  extreme  left  in  figure  5 respresents  the  case  where  bubble  velocity 
is  negligible,  as  with  stoichiometric  propane-air  combustion  in  the  earth's  gravi- 
tational field.  The  data  near  the  left  axis  of  figure  3 represent  this  condition.  If 
a strong  centrifugal  force  field  is  added,  however,  the  bubble  rsces  ahead  of  the 
advancing  turbulent  flamefront,  as  shown  in  the  middle  sketch  of  figure  5.  To  an 
observer  measuring  flame  propagation,  the  rate  would  depend  only  on  centrifugal 
acceleration  intensity.  This  is  the  situation  represented  by  the  data  above  about 
500  g in  figure  3.  The  sketch  on  the  extreme  right  in  figure  5 represents  a case 
where  a centrifugal  force  field  is  applied,  but  the  turbulent  flamespeed  still  ex- 
ceeds the  bubble  velocity.  The  stoichiometric  hydrogen-air  data  of  figure  4 
represent  this  case.  Thus,  it  is  apparent  that  the  higher  of  the  two  velocities, 
either  turbulent  flamespeed  or  bubble  velocity,  determines  the  rate  of  flame 
propagation  through  a combustible  mixture,  and  the  other  has  no  effect. 

It  should  be  noted  that  the  observed  flamespeed  shown  in  figure  4 is  not  the 
true  flamespeed  through  the  cold  mixture,  but  is  the  sum  of  the  true  flamespeed, 
SB,  plus  the  velocity  of  the  unbumed  mixture  through  the  pipe  caused  by  the  ex- 
pansion of  the  burned  gases.  Based  on  the  assumptions  of  no  heat  loss  and  100$ 
combustion  efficiency,  the  expansion  of  the  burned  gas  can  be  calculated  and 
subtracted  from  the  observed  velocity  to  obtain  the  true  velocity  of  the  hot  flame 
bubble  through  the  cold  fuel-air  mixture.  Figure  6 presents  the  bubble  velocity 
calculated  from  the  observed  flame  propagation  rates. 

As  noted  on  the  curve,  above  about  500  g the  bubble  velocity,  SB,  can  be 
represented  by  the  equation 


Sr,  (1.25)  (g  ) 

B 'Kat  igniter' 


Eq  (2) 


for  the  conditions  tested. 
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Figure  6.  Flame  Bubble  Velocity  as  a Function  of  FI)  74018 

Centrifugal  Acceleration 

Below  about  200  g turbulent  flnmespeed  con*  *ols  flame  propagation,  and  the 
rate  is  independent  of  centrifugal  acceleration.  It  is  interesting  to  note  that,  if 
Equation  2 is  extended  down  to  1 g conditions,  the  bubble  velocity  is  1.25  feet 
per  second.  This  velocity  agrees  very  closely  with  published  measurements  of 
air  bubbles  in  water  (Reference  1)  and  flame  hubbies  in  cold  fuel-air  mixtures 
(Reference  2)  near  the  lean  fuel-air  ratio  limit,  where  the  normal  flnmespeed 
is  low  enough  for  buoyant  forces  to  become  dominant  even  at  1 g conditions. 
Additional  calculations  of  true  flnmespeed  based  on  the  measured  rate  of  pres- 
sure rise  near  the  beginning  of  combustion  gave  similar,  but  generally  higher, 
flamespeed  values. 
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Additional  tests  were  run  to  determine  the  pressure  effect  on  flamespeed. 

These  tests,  presented  in  figure  7,  showed  that  the  flamespreading  process  is 
independent  of  pressure  between  75  psia  and  5 psia.  Below  5 psia,  the  experi- 
mental data  deviate  significantly  from  those  at  higher  pressure,  probably  as  a 
result  of  the  increasing  ratio  of  heat  lost  to  heat  generated  as  the  pressure  is 
reduced. 

All  of  the  flamespreading  experimental  data  obtained  under  the  current  con- 
tract are  tabulated  in  Appendix  A. 

a.  Development  of  Analytical  Model 

An  analytical  model  based  on  basic  bubble  mechanics  was  formulated  to 
predict  the  observed  flamespreading  rates.  From  Newton's  law  of  motion  on 
the  flame  bubble,  the  buoyant  force  and  drag  force  must  equate  to  the  mass  times 
the  acceleration  of  the  bubble.  In  equation  form  this  can  be  written  as: 

('b)N)  ($)  * (-r“2)M  (vb)  (*-£)♦  (ab)  (cd)(1/2  ».)  (f)2  E»  <3> 


where 


Density  of  bubble 


Volume  of  bubble 


R Distance  from  center  of  rotation,  + radially  outward 
t - Time  , 

w = Rotation  rate 
fia  - Density  of  cold,  unburned  gas 
Ap  = Projected  area  of  bubble 

- Drag  coefficient  based  on  projected  area. 

Rearranging  terms,  equation  3 becomes: 


Eq  (4) 
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If  it  is  assumed  that  the  flame  rapidly  accelerates  to  a terminal  velocity  where 
buoyant  and  drag  forces  are  equal,  d2R/dt2  becomes  zero  and  the  equation  can 
be  expressed  as: 


Eq  (5) 


Representing  volume,  Vg,  divided  by  the  projected  area,  Ag,  by  a characteristic 
length,  L,  and  the  centrifugal  acceleration  Rw2  by  a,  the  equation  further  refines 


to: 


Eq  (6) 


In  a continuously  flowing  combustion  system,  however,  such  as  a turbojet 
engine  afterburner  and  in  the  combustion  centrifuge,  the  bubble  can  not  accelerate 
to  terminal  velocity  within  the  confines  of  the  combustion  system  and  the  assump- 
tion of  a constant  terminal  velocity  is  not  valid.  In  a flowing  system,  if  the 
centrifugal  force  field  is  generated  by  swirling  the  air  tangentially  about  the 
burner  axis,  the  centrifugal  acceleration  varies  with  radius.  In  the  combustion 
centrifuge,  the  centrifugal  acceleration  also  varies  with  radius.  Because  of  the 
difficulty  in  finding  a closed  solution  to  the  differential  equation,  two,  step- 
integration  computer  programs  were  formulated.  One  program  represents  a 
static  system,  such  as  the  combustion  centrifuge,  and  the  second  represents  a 
flowing  system,  such  as  a turbojet  augmentor.  A description  of  the  programs, 
sample  results,  flow  diagrams,  and  computer  listings  is  presented  in  Appendixes 
B and  C. 

In  equation  4 all  of  the  terms  can  be  explicitly  defined  for  a known  gravita- 
tional field  except  the  term  2V-./A-. C_..  A literature  search  was  made  to  assist 
in  defining  this  term.  Data  in  the  literature  for  air  bubbles  in  liquid  (Refer- 
ence 1)  and  flame  bubbles  in  a cold  fuel-air  mixture  (Reference  2)  show  the  same 
characteristic  shape.  The  bubbles  have  a characteristically  lenticular  shape, 
with  a spherical  leading  surface  and  a flat  or  slightly  concave  back  surface.  The 
shape  is  almost  perfectly  represented  by  a segment  of  a sphere,  as  shown  in  fig- 
ure 8.  In  Reference  3 it  was  shown  that  air  bubbles  form  this  shape  at  Reynolds 
numbers  over  4000,  based  on  bubble  diameter.  A constant  value  of  is  also 
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associated  with  these  Reynolds  numbers.  This  implies  that  a constant  bubble 

half  angle,  0^^,  (figure  S)  is  maintained  for  bubbles  with  Reynolds  numbers 
greater  than  4000. 


Figure  8.  Details  of  Air  or  Flame  Bubble 


FD  89698 


The  drag  coefficient  can  be  calculated  by  equating  the  drag  force  to  the 
buoyant  force  on  the  bubble, 


■ f1 

D \2 


.,2  . 2\ 
v V A 1 


VT^g 


where 


Cq  = Drag  coefficient 
° - Density  of  water 

V = Velocity  of  bubble 
Vg  = Volume  of  l)ubble 
g - Gravitational  acceleration 
A = Radius  of  bubble. 

Rearranging  terms, 

C - VB  g 
D 0.5  V2  a A2 
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The  velocity  of  the  buoble  can  be  calculated  by  using  the  theoretical  model  in 
Reference  1. 

2 / \I/2 

V = l(Ag/sin»maxj  Eq(9) 

This  equation  was  derived  by  holding  the  static  pressure  in  Bernoulli's  equation 
constant  along  the  bubble  surface.  Experimental  data  have  verified  the  model. 

From  geometric  considerations,  the  volume  of  the  bubble  is 

VB  - If  0058  "max  - 008  "max  ' 1 1 • Ec|  <10) 

sm  "max  L J 

Substituting  velocity  and  volume  into  equation  8,  it  becomes 


r _ 18 
D 4 


— COS  0 - COS  o 

s ''max 


max  + 3 


2 

sin  0 


Eq  (11) 


max 


Taking  o = 46°,  C = 0.728.  The  term  2VT}/AT3Cri  can  now  be  calculated, 
max  u a jL5  u 


2V 


B 


2nA * 


3CQs3gmax  " C0S(,max  ' Sj 


.*1 


ABCD  Sin3 1) 


max 


n A2  C 


Eq  (12) 


D 


2V 


B 


abcd 


1.618 


Eq  (13) 


The  term  2V is  a function  of  bubble  radius  A in  an  unconfined  sys- 
tem. In  a confined  system,  such  as  the  combustion  centrifuge,  the  same  rela- 
tionship holds,  except  the  drag  coefficient  is  larger  because  of  pipe  wall  effects 
and  therefore  the  constant  (1.618)  has  a different  value.  From  figure  9,  which  is 
a duplicate  of  figure  10  in  Reference  4,  the  drag  coefficient  in  the  pipe  can  be 

roughly  calculated  to  be  1.  24  for  a 1-inch  radius  bubble.  The  term  2Vt,/A_CT, 

13  13  u 

equals  0.  0303  for  this  drag  coefficient.  Using  this  value  for  2VB/AgCD  in  the 
flame  bubble  model  Jiscussed  in  Appendix  B gave  very  good  agreement  with  ex- 
perimental kua  from  the  combustion  centrifuge.  Figure  10  compares  the 
theoretical  and  experimental  results. 
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Figure  10.  Theoretical  vs  Experimental  Results 
in  Combustion  Centrifuge 
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When  reviewing  the  literature,  it  became  apparent  that  there  were  no  data 
available  on  air  bubbles  at  very  large  Reynolds  numbers.  It  is  hypothesized  that 
a limit  Reynolds  number  is  reached  in  which  either  the  scale  or  intensity  of 
turbulence  becomes  great  enough  to  break  the  bubble  apart.  It  was  reported  in 
Reference  5 that  turbulent  wakes  behind  air  bubbles  do  occur.  Therefore,  the 
beginning  of  turbulence  itself  does  not  break  the  bubble  and  a higher  Reynolds 
number  is  reached  before  bubble  growth  is  limited  by  turbulence. 

An  experimental  investigation  was  undertaken  to  determine  the  limiting 
Reynolds  number.  Single  air  bubbles  of  various  sizes  were  released  in  sea 
water  from  a depth  of  60  feet  and  photographed  as  they  rose  to  the  surface.  (See 
figure  11.)  Air  bubbles  enlarge  as  they  rise  because  of  the  decrease  in  pres- 
sure, and  their  velocity  (and  Reynolds  number)  increase  in  proportion  to  their 
size.  A bubble  breakup  always  occurred  when  a limiting  Reynolds  number  was 
reached.  The  experimental  terminal  velocity  vs  bubble  Reynolds  number  is 
compared  with  theory  in  figure  12.  The  experimental  data  agree  reasonably  well 
with  theoretical  values. 
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Reynolds  Number  Based  on  Bubble  Diameter 

Figure  12.  Terminal  Velocity  vs  Reynolds  FD  89701 

Number  for  Air  Bubbles  Rising 
in  Water 

The  limiting  Reynolds  number  was  found  to  be  equal  to  70,000,  based  on 
bubble  diameter.  In  relatively  unconfined  systems,  such  as  large  swirl  combus- 
tion systems,  the  limiting  Reynolds  number  is  hypothesized  to  define  the  maximum 
bubble  size  if  the  Reynolds  number  of  the  bubble  equals  or  exceeds  70, 000.  Other- 
wise, the  bubble  size  depends  on  the  pilot  dimension. 

2.  Extinction  Limits 

The  second  area  of  investigation  is  the  determination  of  the  centrifugal 
force  limit  at  which  the  flame  is  extinguished.  Because  the  combustion  centrifuge 
is  limited  by  structural  considerations  to  speeds  producing  less  than  5000  g,  it 
was  not  possible  directly  to  measure  the  limiting  g value  over  most  of  the  fuel- 
air  ratio  range.  However,  at  very  lean  and  very  rich  fuel-air  ratios,  the  limiting 
centrifugal  force  values  fell  within  the  centrifuge  operating  range.  At  these  con- 
ditions the  limit  was  extremely  sensitive  to  small  variations  in  fuel-air  ratio  and 
temperature,  and,  as  a result,  the  data  are  somewhat  scattered.  Tests  from  a 
t single  mixture  of  luel  and  air  generally  gave  very  consistent  results,  but  tests 

using  another  batch  mixed  to  the  same  specifications  (*2^)  often  produced  a 
parallel  but  displaced  line  of  data.  Typical  results  for  nominal  0.  0375  and  0.  045 
fuel-air  ratios  are  presented  in  figure  13.  As  shown  in  the  figure,  there  appear 
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to  be  three  different  extinction  mechanisms  controlling  different  regimes  of 
operation.  At  the  bottoms  of  the  two  lines  of  data,  there  is  a low  pressure  limit 
below  which  flame  will  not  propagate  regardless  of  the  centrifugal  force.  This 
pressure  limit  is  the  same  with  the  centrifuge  stationary  and  with  it  operating  to 
produce  the  g values  indicated  on  the  figures.  It  is  believed  thai  this  limit  is 
established  by  flame-quenching  in  the  holes  of  the  turbulence  generator  and  is 
not  pertinent  to  this  investigation.  Above  this  low  pressure  limit  and  below 
about  one  atmosphere,  a '.acond  mechanism  controls  the  flame  extinction  limits. 

In  this  region,  increasing  pressure  permits  operation  to  increasingly  higher 
centrifugal  force  values.  It  is  postulated  that  this  is  a regime  where  viscosity 
effects  are  sufficiently  strong  relative  to  density  differences  so  that  true  flame 
bubbles  controlled  by  buoyant  forces  do  not  become  organized  before  the  turb- 
ulence created  by  their  organization  cools  and  quenches  the  flame.  In  the  major 
regime  of  operation,  above  one  atmosphere  pressure,  the  experimental  data  are 
relatively  consistent  and  agree  quite  well  with  the  results  predicted  ly  the  fol- 
lowing analytical  model. 

a.  Analytical  Extinction  Model 

Figure  14  represents  a flame  bubble  rising  through  a cool  mixture  of  fuel 

and  air  under  the  influence  of  buoyant  forces.  The  heat  generated  in  the  bubble 

can  lie  expressed  as  the  heat  of  combustion  (AH)  multiplied  by  the  mass  of  fuel- 

air  mixture  burned.  This,  in  turn,  can  be  expressed  as  the  density  (p  ) multiplied 

n 

by  the  laminar  flamespeed  (Su)  and  the  instantaneous  area  of  the  bubble  (A). 


Q.n  = (AH)  (Pjj)  (Su)  (A) 


Eq  (14) 


«'«d  aunssaud 


Figure  13.  Lim 
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Figure  14.  Schematic  Representation  of  Flame 
Bubble  Heat  Generation  and  Loss 
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The  heat  lost  from  the  bubble,  neglecting  radiation,  can  be  approximated  by 
the  classical  equation  for  convective  heat  transfer. 


Qout  - (h)  (A)  (AT) 


Eq  (15) 


where  h is  the  convective  heat  transfer  coefficient,  A is  the  heat  transfer  surface 
area,  and  AT  is  the  temperature  difference  between  the  bubble  flame  temperature 
and  the  surrounding  fuel-air  mixture  temperature.  The  convective  heat  transfer 
coefficient,  in  turn,  can  be  evaluated  from 


Nusselt  No. 


(Reynolds  No, ) (Prandtl  No. } 
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where  n is  viscosity,  and  Cg  is  a constant  that  includes  nndprovides  consistent 
units. 

Solving  this  equation  using  physical  properties  based  on  the  cold  fuel-air 
mixture  and  selecting  I)  as  90'/<  oi  the  centrifuge  pipe  diamelet  (based  on  photo- 
graphs in  Reference  1)  results  in  the  predicted  extinction  limits  presented  in 
figure  15.  Very  good  agreement  with  the  experimental  data  is  evident.  It  is 
interesting  to  note  that,  in  addition  to  conversion  factors  to  provide  compatible 
units,  the  constant  Cg,  4. 15  x 10^  in  equation  20,  also  contains  constant  C] 

(from  equation  16).  The  constant  needed  to  obtain  this  agreement  of  experi- 
mental data  and  analytical  predictions  is  0.  0535,  compared  to  the  conventional 


l 
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0. 023  used  for  heat  transfer  between  a fluid  and  a solid  pipe  surface. 
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Because  of  the  ability  to  predict  flame  extinction  from  a cooling  analysis, 
changing  only  the  empirical  constant,  it  appears  that  this  approach  may  also  have 
merit  in  predicting  blowoff  from  conventional  flameholders. 

3.  Reynolds  Number  Effects 

Another  area  of  investigation  involves  the  effects  of  Reynolds  No.  on 
turbulent  flamespeed  in  gravitational  fields  low  enough  so  that  buoyancy  has  no 
effect  on  flame  propagation  rates.  The  rate  of  flame  propagation  and  flame- 
front  shape  were  measured  with  ionization  probes  while  hydrogen-air  or  propane- 
air  mixtures  were  burned  in  polyvinyl  chloride  pipes  of  three  different  diameters. 
The  apparatus  is  shown  in  figure  16.  Each  tube  contained  a sparkplug  located 
behind  a turbulence  generator  to  ensure  the  rapid  generation  of  a turbulent  flame. 
Two  ionization  probes  detected  the  passage  of  the  center  of  the  flamefront  at  a 
point  12  inches  downstream  of  the  turbulence  generator,  and  a second  set  of 
probes  detected  the  passage  of  the  flame  at  five  radial  positions  in  a plane  12 
inches  downstream  of  the  first  set.  Thus,  the  two  sets  of  probes  provided 
measurements  of  both  the  average  speed  between  the  two  sets  and  the  shape  of 
the  flamefront  at  the  second  set. 

To  conduct  a test,  a plate  and  a paper  disk  were  placed  over  the  open  end 
of  the  PVC  tube  (which  was  coated  with  vacuum  grease),  and  a vacuum  was  pulled 
to  remove  all  air,  moisture,  and  products  of  prior  combustion  from  the  rig.  The 
previously  prepared  fuel-air  mixture  was  then  admitted  to  the  pipe  from  a separate 
premixed  supply  tank  until  atmospheric  pressure  was  reached.  The  metal  disk 
was  then  slipped  off,  and  the  paper  disk  was  left  adhering  to  the  vacuum  grease. 

The  sparkplug  was  then  fired,  and  the  flame  passage  was  recorded  at  essentially 
constant  atmospheric  pressure. 

The  observed  flamespeed  was  the  sum  of  the  turbulent  flamespeed  plus  the 
velocity  of  the  cold  mixture  due  to  the  expansion  caused  by  combustion  in  the 
closed  end  of  the  pipe.  If  the  flamefront  is  a flat  plane  and  there  is  no  inefficiency 
or  heat  loss,  the  turbulent  flamespeed  can  be  approximated  from  the  expression, 

ST  T^/t”  Eq  (21) 

where  ST  is  turbulent  flamespeed,  U0  is  the  observed  velocity,  and  T1  and  T2  are 
temperatures  before  and  after  combustion.  This  equation  also  assumes  the  number 
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of  moles  before  and  after  combustion  are  nearly  the  same.  This  is  correct  within 
4$  for  stoichiometric  or  leaner  propane-air  mixtures,  but  a correction  for  the 
almost  15$  change  in  number  of  moles  was  required  for  the  hydrogen-air  mixtures. 


Figure  16.  Tubes  Used  To  Measure  Reynolds  FE  336140 

Number  Influence  on  Flamespeed 
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If  the  flamefront  is  not  a flat  plane,  equation  21  can  be  corrected  for  the 
measured  flamefront  area  as  shown  below: 


Eq  (22) 


where  A^  is  the  flamefront  area  and  A 2 is  the  pipe  cross-sectional  area. 

In  most  of  the  experiments  conducted  A^/Ag  was  so  close  to  1. 0 that  it 
could  be  neglected.  Typical  flamefront  profiles  for  hydrogen-air  and  propane- 
air  mixture.^  are  presented  in  figure  17. 

Using  equation  22  to  calculate  the  turbulent  fk’nespeed,  taking  the  laminar 
flamespeed  from  the  literature,  and  calculating  Reynolds  No.  from  the  cold  gas 
velocity  (UQ  - ST)  resulted  in  the  data  shown  in  figure  18.  The  turbulent-to- 
laminar  flamespeed  ratio  increases  rapidly  as  turbulence  increases  until  a maxi- 
mum is  reached.  Beyond  that  point  the  flame  is  stretched  so  thin  that  some  of  it 
is  cooled  to  extinction,  and  further  increase  in  turbulence  results  in  a decrease 
in  turbulent  flamespeed.  It  is  significant  that  both  propane-air  mixtures  with  a 
stoichiometric  laminar  flamespeed  of  1.3  feet  per  second  and  hydrogen-air 
mixtures  with  a stoichiometric  laminar  flamespeed  of  6.  2 feet  per  second  fit  the 


same  curve. 


PROPANE  - AIR,  F/A  = 0.064 


4.  Application  of  Itesults 


The  end  objective  of  all  research  work  is  the  practical  application  of  the 
researcn  results.  One  practical  application  of  the  increased  flamespeed  re- 
sulting from  centrifugal  effects  is  the  swirl  augmentor.  Augmentors  are  cur- 
rently used  in  aircraft  turbojet  and  turbofan  engines  to  provide  thrust  augmen- 
tation. Conventional  augmentors  use  numerous  flameholdeis,  which  act  as 
ignition  sources  to  burn  the  flowing  fuel-air  mixture  (figure  19).  Large  pres- 
sure losses  (4(/()  and  long  augmentor  lengths  (L/D  ~ 2.5)  are  associated  with 
conventional  augmentors.  As  an  alternative,  the  swirl  augmentor  has  been 
proposed.  Swirl  augmentors  use  precombustion  swirl  vanes  to  swirl  the  flow 
and  produce  a centritugal  field  that  increases  the  burning  rate  of  the  combusti- 
ble mixture.  In  figure  20  a schematic  of  a swirl  augmentor  tested  under  another 
program  is  presented  to  define  the  essential  details.  Fuel  is  inserted  into  the 
swirling  flow  by  sprayrings,  and  the  combustible  mixture  is  ignited  by  a cir- 
cumferential pilot  on  the  outer  wall  of  the  combustion  duct.  Initiating  at  the 
pilot,  the  flame  bubbles  move  quickly  to  the  center  due  to  the  buoyant  effect  of 
the  hot  burned  gases  in  the  cold  unburned  gases,  igniting  the  fuel-air  mixture 
as  they  go.  Since  the  bubble  speed  can  be  many  limes  greater  than  turbulent 
flamespeed  in  conventional  augmentors,  shorter  lengths  or  higher  ourning 
efficiencies  result.  Pressure  losses  are  less  than  with  conventional  augmentors, 
especially  if  movable  vanes  are  used  which  are  set  to  zero  degrees  during  non- 
augmemed  flight.  In  addition,  the  swirl  augmentor  has  the  added  performance 
advantage  of  being  unaffected  by  pressure  and  Mach  number  changes. 


To  permit  practical  application  of  these  principles,  a swirl  combustor 
design  system  has  been  developed,  evaluated,  and  applied.  Details  of  the  design 
system  and  a computer  listing  arc  presented  in  Appendix  C. 


The  model's  prediction  of  radial  flame  movement  is  based  on  equation  4 
and  radial  expansion  of  the  burned  gases.  Axial  movement  is  calculated  from  the 
initial  axial  gas  velocity  increased  by  the  expansion  resulting  from  combustion. 
The  result  is  a flamefront  profile  (similar  to  that  in  figure  25;,  which  is  a func- 
tion of  design  variables. 


Figure  19.  Schematic  of  a Conventional  FD  89702 

Augmentor 


Figure  20.  Schematic  of  a Swirl  Augmentor 


FD  89703 


The  model  was  used  to  predict  the  combustion  efficiency  as  a function  of 
length  for  a 15-incli  diameter  swirl  augmentor  rig,  which  was  tested  under 
another  program.  Unfortunately,  the  maximum  bubble  Reynolds  number  was  not 
reached  in  this  small  rig.  In  order  to  predict  flamespreading  rates,  a bubble 
diameter  was  selected  equal  to  the  pilot  dimension.  A comparison  of  the  pre- 
dicted and  experimental  results  is  presented  in  figure  21.  After  this  verification 
of  the  model,  it  was  used  to  design  a swirl  augmentor  for  an  advanced  version  of 
a large  Air  Force  turbofan  engine.  This  design  is  presented  in  figure  22.  The 
reduced  pressure  loss  and  reduced  sensitivity  of  the  swirl  augmentor  to  augmentor 
pressure  and  gas  velocity  are  predicted  to  produce  a 2 % reduction  in  fuel  con- 
sumption and  to  increase  stability  at  high-altitude,  low-Mach-number  flight  con- 
ditions. For  this  application,  the  augmentor  length  was  not  changed  from  that  of 
the  current  engine,  and  the  benefit  was  taken  as  reduced  pressure  loss  and  fuel 
consumption.  Alternatively,  the  pressure  loss  could  have  been  held  constant  and 
the  benefits  of  reduced  length  and  weight  would  have  resulted,  as  shown  in  fig- 
ure 23. 


Figure  21.  Comparison  of  Experimental  and  FD  89704 

Theoretical  Results  at  Augmentor 
Equivalence  Ratios  Between  0. 95 
and  1.  05 
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Another  application  in  which  the  model  was  used  was  in  sizing  a ramburner 
system  for  an  advanced  high-Mach  aircraft  demonstrator  propulsion  system. 

Since  size  and  weight  were  at  a premium  for  this  application,  every  effort  was 
made  to  minimize  both  of  those  parameters.  As  shown  in  figure  24,  the  swirl 
ramburner  was  designed  using  the  model  and  is  significantly  smaller  in  size  and 
also  weight  when  compared  to  the  more  conventional  preburner-ramburner.  Both 
of  the  ramburner  systems  were  sized  for  the  same  thrust. 


120.5  in. 


Separate 
Preburner  • 
Ramburner 


Figure  24.  Swirl  Burner  Concept  Reduces  Ramburner  FD  90740 
Length  - Ramburners  Sized  for  the  Same 
Thrust 


36 


SECTION  m 
CONCLUSIONS 

A bubble  theory  of  flame  propagation  has  been  developed  and  can  be  suc- 
cessfully used  to  design  swirl  flow  combustion  systems. 

A theory  has  been  developed  that  predicts  with  reasonable  accuracy  the 
extinction  limits  of  a flame  in  a high  centrifugal  force  field. 

The  turbulent  flamespeed  of  a fuel-air  mixture  increases  with  increasing 
Reynolds  number,  until  a maximum  is  reached  beyond  which  flame  stretch 
quenches  flame  ligaments  and  results  in  lower  turbulent  flamespeeds. 


APPENDIX  A 
EXPERIMENTAL  DATA 


Table  A-l.  Centrifuge  Test  Conditions  for  Propane  and  Air 


Run 

f/a 

Pressure, 

psia 

Temperature, 

°F 

rpm 

1 

0.  064 

15 

Unrecorded 

0 

2 

0.  064 

15 

Unrecorded 

0 

3 

0.  064 

15 

Unrecorded 

850 

4 

0.  064 

15 

Unrecorded 

0 

5 

0.  064 

15 

Unrecorded 

0 

6 

0.  064 

15 

Unrecorded 

0 

7 

0.  064 

15 

Unrecorded 

850 

8 

0.  064 

15 

Unrecorded 

850 

9 

0.  064 

15 

Unrecorded 

2,015 

10 

0.  064 

15 

Unrecorded 

1,722 

11 

0.  064 

15 

Unrecorded 

1,570 

12 

0.  064 

15 

Unrecorded 

1,404 

13 

0.  064 

15 

Unrecorded 

1,212 

14 

0.064 

15 

Unrecorded 

990 

15 

0.  064 

15 

Unrecorded 

625 

16 

0.  089 

15 

Uni.  corded 

0 

17 

0.089 

15 

Unrecorded 

993 

18 

0.  089 

15 

83 

1,409 

19 

0.  089 

15 

83 

1,567 

20 

0.089 

15 

98 

0 

21 

0.  089 

15 

99 

1,710 

22 

0.  089 

15 

Unrecorded 

1,820 

23 

0.089 

15 

103 

316 

24 

0.089 

30 

103 

0 

25 

0.  089 

30 

100 

627 

26 

0.  089 

30 

98 

990 

27 

0.  089 

30 

98 

1,813 

28 

0.  089 

5 

98 

1,797 

29 

0.  089 

2 

98 

1,790 

30 

0.  089 

1 

103 

1,491 

31 

0.  089 

1 

97 

0 

32 

0.089 

1 

98 

402 

33 

0.  089 

1 

100 

451 

34 

0.089 

1 

100 

0 

35 

0.089 

1 

102 

0 

36 

0.089 

1 

104 

504 

37 

0.089 

1 

104 

0 

38 

0.089 

2 

104 

1,000 

39 

0.  089 

2 

104 

1,507 

40 

0.  089 

2 

103 

1,602 

11 

0.  089 

2 

106 

1,768 

42 

0.089 

2 

106 

1,765 

13 

0.  089 

1.5 

103 

1,  402 

44 

0.  089 

1.5 

104 

1,369 

45 

0.089 

1.5 

104 

1,298 

Table  A-l.  Centrifuge  Test  Conditions  for  Propane  and  Air  (Continued) 


Run 

f/a 

Pressure, 

Tempera  tur 

psia 

°F 

46 

0.089 

1.5 

101 

47 

0.089 

1.5 

102 

48 

0.089 

1.5 

104 

49 

0.089 

1.5 

105 

50 

0.089 

1.5 

106 

51 

0.089 

1.5 

106 

100 

0.  064 

2.5 

Unrecorded 

101 

0.064 

2.5 

Unrecorded 

102 

0.  064 

5.0 

Unrecorded 

103 

0.064 

2.5 

Unrecorded 

104 

0.064 

5.0 

Unrecorded 

105 

0.  064 

5.  0 

Unrecorded 

106 

0.064 

5.0 

Unrecorded 

107 

0.  064 

5.0 

Unrecorded 

108 

0.  064 

5.0 

Unrecorded 

109 

0.064 

5.0 

Unrecorded 

110 

0.064 

5.0 

Unrecorded 

111 

0.  064 

2.5 

Unrecorded 

112 

0.064 

2.5 

Unrecorded 

113 

0.064 

2.5 

Unrecorded 

114 

0.064 

2.5 

Unrecorded 

115 

0.064 

2.5 

Unrecorded 

116 

0.089 

2 

Unrecorded 

117 

0.089 

2 

Unrecorded 

118 

0.  089 

2 

Unrecorded 

119 

0.089 

2 

Unrecorded 

120 

0.089 

2 

Unrecorded 

121 

0.089 

2 

Unrecorded 

122 

0.089 

2 

Unrecorded 

123 

0.089 

2 

Unrecorded 

124 

0.  099 

5 

Unrecorded 

125 

0.099 

5 

Unrecorded 

126 

0.099 

5 

105 

127 

0.099 

<; 

1' 

103 

128 

0.099 

0 

107 

129 

0.  099 

5 

107 

130 

0.099 

5 

103 

131 

0.099 

5 

Unrecorded 

132 

0.  099 

5 

105 

133 

0.099 

5 

Unrecorded 

134 

0.099 

30 

100 

135 

0.099 

5 

Unrecorded 

136 

0.  099 

5 

Unrecorded 

137 

0.099 

10 

Unrecorded 

138 

0.099 

10 

Unrecorded 

139 

0.099 

10 

Unrecorded 

140 

0.099 

15 

Unrecorded 

141 

0.099 

15 

Unrecorded 

rpm 

0 

490 

746 

876 

1,000 

1,102 

0 

0 

0 

1,200 

0 

101 

625 

1,210 

1,565 

990 

1,714 

625 

1.209 
1,706 

100 

1,571 

0 

625 

1.210 

1,800 

2,000 
1,800 
1,700 

0 

C 

625 

1,980 

1,993 

1,880 

0 

1,936 

0 

2,000 

2,000 

0 

2,100 

2,000 

2,100 

2,150 

2,150 

2,100 

2,100 
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Table  A-l.  Centrifuge  Test  Conditions  for  Propane  and  Air  (Continued) 


Run 

f/a 

Pressure, 

psia 

Temperature, 

°F 

rpm 

142 

0.  099 

15 

Unrecorded 

2,250 

143 

0.099 

15 

Unrecorded 

2,308 

144 

0.099 

5 

Unrecorded 

2,250 

145 

0.099 

30 

Unrecorded 

2,100 

146 

0.105 

15 

Unrecorded 

0 

147 

0. 105 

15 

Unrecorded 

1,600 

157 

0. 112 

15 

105 

0 

163 

0. 112 

5 

Unrecorded 

1,200 

252 

0.064 

15 

95 

500 

258 

0.064 

15 

95 

1,048 

259 

0.064 

5 

90 

500 

265 

0.064 

5 

103 

500 

266 

0.064 

30 

108 

1,500 

267 

0.064 

30 

106 

2,100 

270 

0.  064 

20 

108 

1,900 

272 

0.064 

15 

113 

1,100 

273 

0.064 

15 

113 

1,250 

274 

0.  064 

15 

113 

1,300 

275 

0.064 

15 

115 

1,350 

276 

0.064 

10 

115 

865 

277 

0.  064 

10 

103 

900 

278 

0.064 

2.5 

115 

420 

279 

0.064 

2.5 

116 

375 

280 

0.090 

15 

105 

500 

283 

0.090 

15 

115 

750 

285 

0.090 

30 

114 

1,500 

287 

0.  090 

30 

108 

1,750 

288 

0.  090 

2.5 

110 

240 

289 

0.  090 

2.5 

113 

255 

290 

0.090 

2.5 

118 

277 

291 

0.  090 

2.5 

114 

300 

292 

0.090 

2.5 

118 

320 

295 

0.  090 

10 

118 

510 

296 

0.  090 

5 

125 

340 

298 

0.090 

5 

108 

360 

305 

0.  052 

22.  5 

110 

250 

306 

0.  052 

22.5 

110 

275 

307 

0.052 

22.5 

110 

350 

309 

0.  052 

22.5 

110 

400 

310 

0.  052 

22.5 

110 

450 

315 

0.052 

10 

110 

260 

316 

0.052 

10 

115 

312 

317 

0.052 

10 

115 

360 

319 

0.  052 

10 

115 

400 

324 

0.052 

5 

110 

360 

326 

0.052 

2.5 

115 

306 

329 

0.  058 

10 

110 

550 

336 

0.  058 

2.5 

95 

350 

40 


Table  A-l 


Run 

338 

362 

363 

364 
374 

376 

377 
382 
384 
387 

389 

390 

391 

392 
395 

399 

400 

401 

403 

404 

405 

406 

407 

408 


. Centrifuge  Test  Conditions  for  Propane  and  Air  (Continued) 


Pressure, 


f/a 

psia 

0.  061 

5 

0.  064 

10 

0.  064 

10 

0.  064 

10 

0.  064 

20 

0.  064 

20 

0.  101 

15 

0. 101 

15 

0. 101 

15 

0.  101 

30 

0.  101 

5 

0.  101 

4.8 

0.  101 

4.  8 

0.  101 

4.8 

0.  101 

4.7 

0.  101 

20 

0.  101 

20 

0.  101 

20 

0.  101 

30 

0.  101 

30 

0.  101 

30 

0.  101 

30 

0.  101 

28.  3 

0.  101 

20 

Temperature, 


°F 

rpm 

90 

462 

110 

700 

110 

750 

100 

775 

90 

0 

90 

1,200 

95 

225 

95 

195 

85 

195 

85 

310 

85 

121 

85 

125 

85 

150 

85 

200 

100 

190 

110 

290 

115 

315 

110 

327 

115 

439 

115 

460 

115 

505 

120 

550 

120 

550 

120 

330 
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Table  A-2.  Centrifuge 


Run 

Probe 

Ft  From  Center 

Time 

Sec  After  Soark 

Probe 

Ft  From  Center 

Time 

Sec  After  Spark 

1 

2.416 

0.  01676 

2.333 

0.01756 

2 

2.416 

0.01627 

2.333 

0.01711 

4 

2.416 

0.01579 

2.333 

0. 01665 

5 

2.416 

0. 01634 

2.333 

0.01707 

6 

2.416 

0. 01699 

2.333 

0.01765 

7 

2.416 

0. 02005 

2.000 

0. 02518 

9 

2.416 

0.  00468 

2.  000 

0.00678 

10 

2.416 

0.00669 

2.000 

0. 00810 

11 

2.416 

0.00681 

2.000 

0.00831 

12 

2.416 

0.00760 

2.  000 

0. 00922 

13 

2.416 

0. 00805 

2.000 

0. 00994 

14 

2.416 

0.00909 

2.000 

0.00140 

15 

2.416 

0.01138 

2.  000 

0.01449 

16 

2.416 

0.  01423 

2.000 

0.01768 

17 

2.416 

0.  00841 

2.000 

0.01061 

18 

2.416 

0.  00729 

2.000 

0. 00898 

19 

2.416 

0.  00721 

2.  000 

0. 00861 

20 

2.416 

0. 01550 

2.000 

0.01957 

21 

2.416 

0.00630 

2.000 

0. 00769 

22 

2.416 

0.  00652 

2.  000 

0.00786 

23 

2.000 

0. 02055 

1.916 

0.02141 

24 

2.416 

0.  02092 

2.250 

0.  02276 

25 

2.416 

0. 01352 

2.  250 

0.01406 

26 

2.416 

0.00835 

2.  250 

0.  00933 

27 

2.250 

0. 00639 

2.000 

0.00710 

28 

2.416 

0.  00552 

2.  250 

0.00638 

29 

2.416 

0.  00656 

2.  250 

0.  00803 

30 

2.416 

0.01006 

2.250 

0.  01387 

31 

2. 416 

0. 02767 

2.250 

0.  03307 

32 

2.416 

0.01957 

2.  250 

0.  02798 

33 

2.416 

0.01964 

2.250 

0.  02711 

34 

2.416 

0. 08446 

2.  250 

0.11286 

36 

2.416 

0.01756 

2.  250 

0.  02232 

37 

2.416 

0. 03054 

2.250 

0.  03643 

38 

2.416 

0. 010/6 

2.250 

0.01233 

33 

2.416 

0.00914 

2.250 

0.01109 

40 

2.416 

0.00850 

2.250 

0.01383 

4i 

2.416 

0. 00626 

2.250 

0.00790 

42 

2.416 

0.00838 

2.250 

0.01089 

43 

2.416 

0.01095 

2.250 

0.01500 

44 

2.416 

0.01042 

2.250 

0.01349 

45 

2.416 

0.01041 

2.250 

0.01367 

46 

2.416 

0.02108 

2.250 

0.  02401 

47 

2.416 

0.01538 

2.  250 

0.01763 

48 

2.416 

0.01231 

2.250 

0.01428 

49 

2.416 

0.01170 

2.250 

0. 01455 

50 

2.416 

0. 01253 

2.250 

0.  01485 

51 

2.416 

0.01413 

2.250 

0.01850 

100 

2.416 

0.01518 

2.  250 

0.01738 

Probe  Tini 

Ft  I rom  Center  Sec  Aftei 


2.000 
2.  000 
2.  000 

2.000 

2.000 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 

1.916 
1.583 

2.000 

2.000 

2.  000 

1.916 
2.  000 
2.  000 
2.  000 
2.  000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 


0.02 
0.02 
0.01 
0.  02 
0.02 
0.02 
0. 00 
0.00 
0.  00 
0.00 
0.01 
0.  01 
0.  01 
0.01 
0.01 
0.  00 
0.00 
0. 02 
0.00 
0.00 
0.02 
0.02 
0.01 
0.01 
0. 00 
0.00 
0.00 
0.01 
0.  03 
9.03 
0.03 
0. 12 
0.02 
0.04 
0.01 
0.01 
0.  01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.  03 
0.01 
0.011 
o.or 
0.01 
0.  02, 
0.02' 


OM| 


ft 


,1 


[Table  A-2.  Centrifuge  Test  Data  for  Propane  and  Air 


Probe 
tom  Center 


.000 
.000 
. 000 
.000 
.000 
; 91G 
.916 
.916 
,.916 
:.916 
.916 
.816 
916 
.916 
.916 
.916 
,916 
.916 
'.916 

r.  9i6 

.583 

.000 

.000 

.000 

.916 

tOOO 

.000 

.000 

.000 

.000 

;000 

,000 

s000 

.000 

-Ann 

p vvv 
>000 
,000 
,000 
.000 
,000 
,000 
,000 
;ooo 
000 
.000 
000 
,000 
000 
000 


Time 

Sec  After  Spark 

Probe 

Ft  From  Center 

Time 

Sec  After  Spark 

Probe 

Ft  From  Center 

Time 

Sec  After  Spark 

Probe 
Ft  From  C< 

0.02085 

1.916 

0.02185 

1.583 

0.02585 

1.500 

0. 02060 

1.916 

0.02175 

1.583 

0. 02693 

1.500 

0. 01982 

1.  916 

0. 02067 

1.583 

0.02537 

1.500 

0. 02073 

1.916 

0.02171 

1.583 

0. 02878 

1.500 

0.02114 

1.916 

0.02217 

1.583 

0.  02825 

1.  500 

0. 02601 

1.583 

0.02912 

1.500 

0.02971 

0. 00696 

1.583 

0. 00766 

1.500 

0. 00784 

0.00831 

1.583 

0.00908 

1.500 

0. 00826 

0. 00855 

1.583 

0.00934 

1.500 

0. 00952 

0. 00946 

1.583 

0. 01036 

1.500 

0. 01060 

0.01024 

1.583 

0.01128 

1.500 

0.01152 

0.01171 

1.583 

0.01287 

1.500 

0.01317 

0.01491 

1.583 

0.01653 

1.500 

0.  01689 

0.01839 

1.583 

0. 02143 

1.500 

0.02220 

0.01091 

1.583 

0.01213 

1.500 

0.  01244 

0.00922 

1.  583 

0.01006 

1.500 

0.01024 

0.00885 

1.583 

0. 00964 

1.500 

0. 00982 

0.02018 

1.583 

0. 02383 

1.500 

0. 02462 

0.00793 

1.583 

0. 00854 

1.500 

0.  00872 

0. 00804 

1.583 

0. 00865 

1.500 

0.00878 

0. 02436 

1.500 

0. 02497 

0.02528 

1.916 

0.02613 

1.666 

0. 02859 

1.500 

0.01570 

1.916 

0.01618 

1 . 666 

0.01745 

1.500 

0.01030 

1.916 

0.01055 

1. 666 

0.01140 

1.500 

0. 00728 

1.666 

0. 00781 

1.500 

0.00817 

0.00718 

1.916 

0. 00736 

1-666 

0.00798 

1.500 

0.00987 

1.916 

0.01042 

1.666 

0.01171 

1.500 

0.01816 

1.916 

0.01933 

1.666 

0.02301 

1.500 

0.  03969 

1.916 

0.04190 

1.666 

0. 04485 

1.500 

0. 03380 

1.916 

0.03491 

1.066 

0. 03767 

1.500 

0. 03301 

1.916 

0.0341G 

1.666 

0.03723 

1.500 

0. 12946 

1.916 

0. 02598 

1.916 

0. 02695 

1.666 

0. 02970 

1.500 

0. 04458 

1.916 

0.04571 

1.666 

0. 04988 

1.  500 

0.01121 

1. 916 

n 01,171 

1.  666 

0. 01622 

1.500 

0.01317 

1.916 

0.01371 

1.666 

0.01512 

1.500 

0.01383 

1.  916 

0.01455 

1.666 

0.01671 

1.500 

0.01012 

1.916 

0.01080 

1.666 

0.01264 

1.500 

0.01407 

1.916 

0.01491 

1.666 

0.01706 

1.500 

0.01917 

1.916 

0.02042 

1.666 

0. 02417 

1.500 

0.01687 

1.916 

0.01766 

1.666 

0.01963 

1.500 

0. 01675 

1.916 

0.01757 

1.666 

0.01970 

1.500 

0. 03108 

1.916 

0.03479 

1.666 

0. 04527 

0.01982 

1.916 

0.02041 

1.666 

0. 02207 

1.500 

0. 01642 

1.916 

0.01699 

1.666 

0.01855 

1.500 

0.01752 

1.916 

0.01842 

1.666 

0. 02048 

1.500 

0.01779 

1.916 

0.01859 

1.666 

0.  02055 

1.500 

0. 02359 

1.  916 

0.02487 

1.666 

0. 02814 

1.500 

0. 02095 

1.916 

0.02327 

1 .666 

0. 02827 

1.500 

Id  Air 


V 


Time 

Probe 

Time 

■After  Spark 

Ft  From  Center 

Sec  After  Spark 

Ft  From  Center 

0. 02185 

1.583 

0. 02585 

1.  500 

0.02175 

1.583 

0.02693 

1. 500 

0. 02067 

1.583 

0. 02537 

1.500 

0. 02171 

1.583 

0.  02878 

1.500 

0. 02217 

1.583 

0. 02825 

1.500 

0.02912 

1.500 

0.02971 

0.  00766 

1.500 

0. 00784 

0.00908 

1.500 

0.00826 

0.  00934 

1.500 

0.00952 

0.  01036 

1.500 

0.01060 

0.  01128 

1.500 

0.01152 

0.01287 

1.500 

0.01317 

0.01653 

1.500 

0.01689 

0.  02143 

1.500 

0.02220 

0.01213 

1.500 

0.  01244 

0.01006 

1.500 

0.  01024 

0.  00964 

1.500 

0. 00982 

0.  02383 

1.500 

0.02462 

0.  00854 

1.500 

0. 00872 

0. 00865 

1.500 

0.00878 

0. 02497 
0. 02613 

1.  666 

0.  02859 

1.500 

0.01618 

1.  666 

0 01745 

1.500 

0.01055 

1.666 

0.01140 

1.500 

0.00781 

1.  500 

0. 00817 

0. 00736 

1.666 

0.00798 

1.500 

0.01042 

1. 666 

0.01171 

1.500 

0.  01933 

1.666 

0. 02301 

1.500 

0. 04190 

1.666 

0. 04485 

1.500 

0.03491 

1.666 

0. 03767 

1.500 

0. 03416 

1.666 

0.03723 

1.500 

0.  02695 

1. 666 

0.02970 

1.500 

0. 04571 

1.666 

0.04988 

1.500 

0.01471 

1.666 

0.01622 

1.500 

0.  01371 

1.666 

0.01512 

1.500 

0.01455 

1. 666 

0.01671 

1.500 

0.01080 

1.666 

0.01264 

1.500 

0.01491 

1.666 

0.  01706 

1.500 

0. 02042 

1.666 

0. 02417 

1.500 

0.01766 

1.666 

0.01963 

1.500 

0.01757 

1.666 

0.01970 

1.500 

0. 03479 

1.666 

0.  04527 

0. 02041 

1.666 

0.  02207 

1.500 

0.  01699 

1.666 

0.01855 

1.500 

0.01842 

1.666 

0. 02048 

1.500 

0.01859 

1.666 

0.  02055 

1.500 

0.  02487 

1. 666 

0. 02814 

1.500 

0. 02327 

1.666 

0,  02827 

1.500 

Sec 


Time 

After  Spark 

0.02701 
0.02765 
0 ^2610 
0.02927 
0. 02934 


0.03018 

0.01818 

0.01195 

0.00840 
0.01257 
0. 02472 
0.  06086 
0. 03896 
0.  03880 

0.  03122 
0.  06851 
0.01703 
0.  01591 
0.0J766 
0.01350 
0.01838 
0.02601 
0.  02098 
0. 02083 

0.  02302 
0.01942 
0.02152 
0. 02153 
0. 02964 
0. 03137 


42 


Table  A-2.  Centrifuge  Test 


Run 

Probe 

Ft  From  Center 

Time 

Sec  After  Spark 

Probe 
Ft  F'rom  C 

101 

2.416 

0.01726 

2.  250 

102 

2.416 

0.01445 

2. 2"' 

103 

2.250 

0. 01045 

2.  000 

104 

2.250 

0.01560 

2.000 

105 

2.250 

0.01494 

2.  000 

106 

2.250 

0.01283 

2.  000 

107 

2.250 

0. 00928 

2. 166 

108 

2.416 

0. 00744 

2.250 

109 

2.416 

0. 00936 

2.250 

110 

2.416 

0. 00645 

2.250 

111 

2.416 

0.01257 

2.250 

112 

2.416 

0. 00994 

2.250 

113 

2.416 

0.00878 

2.250 

114 

2.416 

0.01620 

2.250 

115 

2.416 

0. 00825 

2.250 

116 

2.416 

0.01329 

2.250 

117 

2.416 

0.01134 

2.25C 

118 

2.416 

0.  00829 

2.250 

119 

2.416 

0.00729 

2.250 

120 

2.416 

0.00726 

2.250 

121 

2.416 

0.00707 

2.  250 

122 

2.416 

0.00866 

2.250 

123 

2.416 

0.01418 

2.250 

124 

2.416 

0.01834 

2.  250 

125 

2.250 

0.01585 

2.000 

126 

2.250 

0.00933 

2.000 

127 

2. 250 

0.01207 

2.000 

128 

2.416 

0.00831 

2.250 

129 

2.416 

0.01935 

2.250 

130 

2.416 

0. 00747 

2.  250 

131 

2.333 

0. 02085 

2.250 

132 

2.333 

0.00988 

2.250 

133 

2.333 

0.01006 

2.250 

134 

2,333 

0. 04560 

2.250 

135 

2.  333 

0.00915 

2.  083 

136 

2.333 

0. 00825 

2.  250 

137 

2,333 

0.00821 

2.  250 

138 

2. 333 

0. 00855 

2.  250 

139 

2.250 

0.00910 

2.  083 

140 

2.250 

0. 00846 

2.  083 

141 

2.083 

0. 00938 

1.666 

142 

2.250 

0. 00841 

2.  083 

143 

2.250 

0. 00785 

2.083 

144 

2.250 

0. 00828 

2.  083 

145 

2.250 

0. 00805 

2.  083 

146 

2.250 

0. 03059 

2.  083 

147 

2.250 

0.01373 

2.  083 

157 

2.250 

0.05018 

2.  083 

163 

2.083 

0.02091 

1.833 

Time 

See  After  Snarly 

0.01940 
0. 01642 
0 01190 
0.01815 
0. 01768 
0. 01446 
0.01000 
0. 00835 
0. 01064 
0.  00733 
0. 01433 
0.01195 
0.01071 
0.01834 
0.01024 
0.01537 
0.01311 
0.01012 
0. 00928 
0. 00939 
0. 00909 
0.01146 
0. 01636 
0. 02123 
0.01805 
0.01218 
0.01549 
0.01028 
0. 02207 
0. 00932 
0. 02220 
0.01176 
0.01140 
0.04765 
0.01230 
0. 0092? 

0. 00911 
0. 00936 
0.01006 
0. 00931 
0.01084 
0.  00935 
0.  00853 
0.  00911 
0.  00858 
0.  03391 
0.  01596 
0. 05710 
0. 02491 


Probe 

Ft  From  Center 


2.000 

2.000 

1.  916 
1.916 
1.916 
1.916 
1.  666 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
1.916 

2.  000 
1.666 
2.000 
2.000 
2.000 
1.916 
1.666 
1.583 
2.000 
2.000 
2.000 
2.083 
2.083 
2.083 
2.083 
1.833 
2.  083 
2.083 
2.083 
1.833 
1.833 
1.583 
1.666 
1.666 
1.833 
1.833 
1.833 
1.833 
1.833 
1.666 


$ 

Sec  Af 


0. 

0. 

0. 

0. 

0., 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

O.i 

O.i 

0.< 

o.< 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.< 

0.1 

0J 

0.< 

0.< 

o.< 

0.( 

o.< 

o.< 

o.< 

o.< 

0.' 

0.( 

o.« 

0.< 

o.< 

0.( 

0.( 

o.< 

0.( 

0.< 

0.( 

0.( 

0.( 

o.c 


ble  A-2.  Centrifuge  Test  Data  for  Propane  and  Air  (Continued) 


Probe 

Time 

Probe 

Ft  From  Center 

Sec  After  Soark 

Ft  From  Center 

2.000 

0.02274 

1.916 

2.000 

0.01896 

1.916 

] . 91G 

0.01230 

1.666 

1.916 

0.01923 

1. 666 

1.916 

0.01890 

1.666 

1.  916 

0.01494 

1.666 

1.666 

0.01151 

1.583 

2.000 

0.00927 

1.916 

2.000 

0.01187 

1.916 

2.  000 

0.00826 

1.916 

2.000 

0.01608 

1.916 

2.000 

0.01385 

1.916 

2.000 

0.01296 

1.916 

2.000 

0. 02110 

1.916 

2.000 

0.01217 

1.916 

2.000 

0.01866 

1.916 

2.000 

0.01500 

1.916 

2.000 

0.01244 

1.916 

1.916 

0. 01253 

1.666 

2.000 

0.01195 

1.666 

1.666 

0. 01482 

2.000 

0.01482 

1.  666 

2.000 

0. 01933 

1.916 

2.000 

0. 02509 

1.916 

1.916 

0.01872 

1.666 

1.666 

0.01521 

1.583 

1.583 

0.01848 

2.000 

0.01298 

1.583 

2.000 

0.  02669 

1.916 

2.000 

0.01179 

1.583 

2.083 

0. 02451 

1.916 

2.  083 

0.01406 

1.  83  ** 

2.083 

0.01297 

1.833 

2.083 

0. 05133 

1.833 

1.833 

0. 01442 

1.666 

9 009 

“ • v w 

A A 1 AAA 

1 ooo 

* ViVUU 

X#  OOO 

2.  083 

0.01042 

1.  833 

2.083 

0.01006 

1.  81-.3 

1.833 

0.01127 

1.66  6 

1.833 

0.01034 

1.666 

1.583 

0.01107 

1.666 

0.01082 

1.583 

1.666 

0.01000 

1.583 

1.833 

0.01041 

1.666 

1.833 

0. 00923 

1.666 

1.833 

0.  03929 

1.583 

1.  833 

0.01753 

1.583 

1.833 

0.06562 

1.666 

1.666 

0.  02754 

1.583 

Time 

Probe 

Time 

Sec  After  Soark 

Ft  From  Center 

Sec  After  Spark 

0.  02405 

1.666 

0. 02881 

0. 02O17 

1.666 

0. 92364 

0.01333 

1.500 

0.01397 

0. 02363 

1.583 

0. 02524 

0.  02341 

1. 583 

0.  02457 

0.01620 

0.01175 

1.583 

0.01663 

0.  00951 

1.666 

0.01018 

0.  01228 

1.666 

0.01327 

0.00855 

1.666 

0. 00913 

0.01661 

1.666 

0.01801 

0. 01444 

1.666 

0.01586 

0.01357 

1.666 

0.01520 

0. 02227 

1.666 

0.  02669 

0.01277 

1.666 

0.01410 

0.  02043 

1.666 

0.  02598 

0.01555 

1.666 

0.01701 

0.01317 

1.666 

0.01488 

0.01434 

1.500 

0.01530 

0.01463 

1.500 

0.01573 

0. 01756 

0.  02061 

1.666 

0.  02545 

0. 02650 

1.666 

0. 03092 

0. 02049 
0. 01582 

1.583 

0. 02104 

0.01579 

0. 02882 
0. 01426 

1.666 

0.  03456 

0. 02744 

1.666 

0.  03220 

0.01659 

1.666 

0.01776 

0. 01506 

1.666 

0.01610 

0. 05572 

1.666 

0. 05910 

0.01564 

1.583 

0.01612 

0. 01295 

1. 66b 

0.01404 

0,01179 

1.666 

0.01238 

0.01122 

1.666 

0.01203 

0.01187 

1.583 

0.01223 

0.01091 

0.01106 

0.01025 

1.583 

0.01114 

0.01112 

1.583 

0.01136 

0.  00953 

1.583 

0.  00964 

0. 04538 

1.500 

0. 04728 

0.01867 

1.500 

0.01910 

0.  07438 
0.  02374 

1.583 

0.  07911 

- 


1 


“ (Continued) 

Time 

ter  Sec  After  Spark 

0. 02405 
0. 02017 
0. 01333 
0. 02363 
0.02341 
0.01620 
0.01175 
0. 00951 
0.01228 
0. 00855 
0.01661 
0.01444 
0.01357 
■ 0.02227 

0.01277 
[ 0.02043 

’ 0. 01555 

0 01317 
0.01434 
0. 01463 

0.01756 
0.02061 
0. 02650 
0. 02049 
0.01582 

0.01579 

f 0. 02882 

0. 01426 
s 0. 02744 

0.01659 
: 0. 01506 

0.05572 

‘ A A1  r/*  4 

' u. uiout 

0.01295 

0.01179 

0.01122 

0.01187 

0.01091 

0.01106 
0.01025 
0.01112 
0. 00953 
0. 04538 
0. 01867 
0. 07438 
0. 02874 


Probe 

Ft  From  Center 

1.666 

1.666 

1.500 

1.583 

1.583 

1.583 

1.666 
1.666 
1.666 
1.666 
1.666 
1.666 
1.666 
1.666 
1.666 
1.666 
1. 666 
1.500 
1.500 


1. 666 
1.666 
1.583 


1.666 

1.666 
1.666 
1.666 
1.666 
l.  563 
1.666 
1.666 
*.666 
1.583 
1.563 


1.583 

1.583 

1.500 

1.500 

1.583 


Time 

Sec  After  Spark 

0. 02881 
0. 02364 
0.01397 
0. 02524 
0. 02457 
0.01663 

0.01018 

0.01327 

0.00913 

0.01801 

0.01586 

0.01520 

0.02669 

0.01410 

0.02598 

0.01701 

0.01488 

0.01530 

0.01573 


0.02545 
0.  03«-92 
0. 02104 


0. 03456 

0. 03220 
0.01776 
0.01610 
0. 05910 
0. 01612 
0. 01404 
0.01238 
0.01203 
0.01223 
0.01114 


0.01136 
0. 00964 
0.04728 
0.01910 
0.07911 


Probe 

Ft  From  Center 

*.500 

1.500 


1.500 

1.500 

1.500 

1.500 

1.500 

1.500 

1.500 

1.500 

1.500 

1.500 

1.500 


1.583 

1.583 


1.583 

1.583 

1.583 


Time 

Sec  After  Spark 

0. 03167 
0.  02607 


0.01055 
0.01386 
0.00953 
0.01889 
0. 01675 
0.01602 
0.02939 
0. 01494 
0. 02835 
0.01793 
0.01591 


0.01829 
0. 01680 


0.01446 

0.01268 

0.01227 
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Table  A-2.  Centrifuge  Test  Data 


. . 


Probe 

Time 

Probe 

Time 

Probe 

T: 

Run 

Ft  From  Center 

Sec  After  Spark 

Ft  From  Center 

See  After  Spark 

Ft  From  Center 

Sec  Afte 

252 

2.250 

0. 1356 

2.  083 

0. 1270 

1.833 

0. 1C 

258 

2.250 

0.2421 

2.  083 

0. 2226 

1.  833 

0.  IE 

259 

2.C83 

0.2319 

1.833 

0. 1361 

1.666 

0. 1( 

265 

2.250 

0.2437 

2.083 

0.2340 

1.833 

0.  IE 

266 

2.250 

0.1275 

2.083 

0.1203 

1.833 

0. 1C 

267 

2.250 

0.0709 

2.  083 

0. 0636 

1.833 

0.  0E 

270 

2.250 

0.1155 

2.  083 

0. 0971 

1.833 

0.  07 

272 

2.250 

0.1139 

2.  083 

0. 1053 

1.  833 

0.  01 

273 

2.250 

0. 1550 

2.  083 

0. 1407 

1.833 

0. 1( 

274 

2.250 

0. 1688 

2.083 

0. 1522 

1.833 

0.  11 

275 

1.833 

0. 1293 

1.666 

0. 1058 

1.583 

0.  Of 

276 

2.250 

0. 1927 

2.083 

0. 1799 

1.833 

0.  IE 

277 

2.250 

0.2248 

2.  083 

0. 2049 

1.833 

0. 16 

278 

1.833 

0. 1710 

1.666 

0. 1097 

1.583 

0.  0* 

279 

2.250 

0. 2447 

2.083 

0. 2096 

1.833 

0.  11 

280 

2.250 

0.2195 

2.  083 

0. 2093 

1.833 

0. If 

283 

2.083 

0.  2757 

1.833 

0. 2296 

1.666 

0.  2( 

285 

2.250 

0.0951 

2.  083 

0. 0907 

1.  833 

0.  07 

287 

2.250 

0.  1340 

2.  083 

0. 1257 

1.  833 

0.  01 

288 

2.250 

0.  1609 

2.083 

0. 1560 

1.  833 

0. 11 

289 

2.250 

0. 1850 

2.  083 

0. 1792 

1.833 

0.  V 

290 

2.250 

0. 1786 

2.083 

0. 1728 

1.  833 

0. 1J 

291 

2.250 

0.1942 

2.083 

0. 1883 

1,833 

0.17 

292 

2.250 

0.2053 

2.083 

0. 2000 

1.833 

0. 17 

295 

2.250 

0. 2396 

2.  083 

0. 2333 

1.  833 

0.  21 

296 

2.250 

0.  2583 

2.083 

0.2490 

1.833 

0.2: 

298 

2.250 

0. 1903 

2.083 

0. 1850 

1.  833 

0.  1( 

305 

2.  083 

0.2434 

1.  833 

0.2078 

1.666 

0. 1( 

306 

2.083 

0. 1874 

1.833 

0. 1587 

1.666 

0. 11 

307 

2.083 

0.2176 

1.833 

0.1878 

1.666 

0. 1( 

309 

2.083 

0.  2255 

1.833 

0. 1966 

1.666 

o.  r 

310 

2.083 

0.2616 

1.  833 

0. 2217 

1.666 

0.  2( 

315 

2.083 

0. 2328 

1.833 

0. 1971 

1.666 

0. 1! 

316 

2.083 

0.2683 

1.  833 

0.2000 

1. 666 

0.  1! 

317 

2.083 

0. 2623 

1.  833 

0. 2034 

1.666 

0. 1( 

319 

2.083 

0.  3024 

1.833 

0. 2327 

1.666 

0. 15 

324 

1.666 

0. 1754 

1.583 

0. 1556 

326 

1.666 

0. 2245 

1.583 

0. 1769 

329 

2.083 

0. 1995 

1.833 

0.  1652 

1.666 

o.  i: 

336 

1,  833 

0.2049 

1.666 

0.1151 

1.583 

o.  i: 

338 

1.666 

0. 1296 

1.583 

0. 1087 

362 

2.083 

0. 1351 

1.  833 

0. 1139 

1.666 

0.  0‘ 

363 

2.083 

0.  2401 

1.  833 

0.2043 

1.666 

0.  1( 

364 

2.083 

0.2153 

1.833 

0.  1756 

1.666 

0. 1' 

374 

2.083 

0.  0748 

1.  833 

0.0592 

1.666 

0.  Of 

376 

2.083 

0. 1700 

1.833 

0.1461 

1.666 

0. 11 

377 

2.083 

0. 1198 

1.833 

0.  0981 

1.666 

0.  01 

382 

2.083 

0. 2757 

1.833 

0. 2583 

1.666 

C 2‘ 

384 

2.083 

0.  3816 

1.833 

0.  3701 

1.666 

0.  3! 

Table  A-2.  Centrifuge  Test  Data  for  Propane  and  Air  (Continued) 


Probe 


Time 


Probe 


Time 


Probe 


Time 


J 


Air  (Continued) 


'I  ime 

Sec  After  Spark 


Probe 

Ft  From  Center 


Time 

Sec  After  Spark 


Probe 

Ft  From  Center 


Time 

Sec  After  Spark 


Table  A -2.  Centrifuge  Test  Data  fo 
Probe  Time  Probe  Time  Probe  Ti 


Run 

Ft  From  Center 

Sec  After  Spark 

Ft  From  Center 

See  After  Spark 

Ft  From  Center 

Sec  Afte 

387 

2.083 

0.5268 

1.833 

0. 5072 

1.666 

0.4 

389 

2.083 

0. 2726 

1.833 

0. 1678 

1.666 

0.1 

390 

2.083 

0.0815 

1.833 

0.  0629 

1.666 

0.  0 

391 

2.083 

0. 1585 

1.  833 

0. 1268 

1.666 

0. 1 

392 

2.083 

0.2178 

1.  833 

0. 1966 

1.666 

0. 1 

395 

2.083 

0. 2362 

1.833 

0. 2193 

1.666 

0.  r 

399 

2.083 

0.3818 

1.833 

0. 3636 

1.666 

o.  3: 

400 

2.083 

0.3101 

1.833 

0. 2859 

1.666 

0.  2' 

101 

2.083 

0. 3302 

1.833 

0. 3142 

1.666 

0. 2! 

403 

2.083 

0. 2353 

1.833 

0. 2092 

1.666 

0. 1! 

404 

2.  083 

0. 2263 

1.833 

0. 1985 

1.666 

0.  li 

405 

2.  083 

0. 2254 

1.833 

0.  2124 

1.666 

0. 2( 

406 

2.083 

0. 2039 

1. 833 

0.1859 

1.666 

0.  It 

407 

2.083 

0.2185 

1.833 

0. 1924 

1.666 

0. 1( 

408 

2.083 

0. 2889 

1.833 

0.  2705 

1.666 

0. 2‘ 

p 


f 


,nd  Air  (Continued) 


Probe 

From  Center 

1. 583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 


I 


Time  Probe  Time  Probe  Time 


Run 

220 
221 
* 222 

223 

224 

225 

226 

227 

228 
229 
231 

235 

236 

240 

241 
243 
245 
2 V' 
248 


Table  A-3.  Centrifuge  Test  Conditions  for  Hydrogen  and  Air 


Pressure, 

Temperature, 

Pipe 

I/a 

psia 

°F 

Size 

0.029 

5 

105 

0 

0.029 

5 

105 

50 

0.029 

5 

110 

100 

0.029 

5 

106 

150 

0.029 

5 

110 

200 

0.029 

5 

103 

300 

0.029 

5 

110 

600 

0.029 

5 

105 

1,200 

0.  029 

5 

112 

2, 112 

0.029 

5 

108 

2,110 

0.  013 

5 

Unrecorded 

0 

0.013 

5 

98 

750 

0.013 

5 

96 

1,000 

0.013 

15 

106 

0 

0.013 

15 

106 

1,500 

0.013 

10 

105 

0 

0.013 

10 

110 

1,500 

0.013 

10 

no 

1,750 

0.013 

10 

100 

1,750 

46 


Run 

220 

221 

222 

223 

224 

225 
22(5 

227 

228 
229 

234 

235 

236 

240 

241 
243 

245 

246 
248 


Table  A-4.  Centrifuge  Test  Data  f< 


Probe 

Ft  From  Center 

2.250 

2.250 

2.250 

2.250 

2.250 

2.250 

2.250 

2.250 

2.250 

2.250 

2.250 

2.  250 

2.250 

2.250 

2.250 

2.250 

2.250 

2.250 

2.250 


Time 

Sec  After  Spark 

U.  01875 
0.01962 
0.01947 
0. 01744 
0.01963 
0.01768 
0.01744 
0. 02305 
0. 04671 

0.  1355 
0. 1109 
0. 1721 
0. 0867 
0.0795 
0. 1026 
0. 0998 
0. 0998 
0. 1056 


Probe 

Ft  From  Center 

2.  083 
2.083 
2.083 
2.  083 
2.083 
2.083 
2.  083 
2.083 
2.  083 
2.  083 
2.  083 
2.  083 
2.083 
2.  083 
2.  083 
2.  083 
2.  083 
2.  083 
2.  083 


Time 

Sec  After  Spark 

0.01635 
0.01722 
0.01731 
0.01549 
0.01720 
0.01561 
0.01561 
0. 02110 
0. 04500 

0. 1254 
0. 0989 
0. 1612 
0. 0801 
0. 0746 
0. 0911 
0.  0954 
0.  0960 
0.  1016 


Probe 

Ft  From  Center 

1.833 
1.833 
1.833 
1.833 
1.833 
1.833 
1.833 
1.833 
1.833 
1.833 
1.833 
1.833 
1.666 
1.833 
1.833 
1.833 
1.  833 
1.833 
1.833 


Itrifuge  Test  Data  for  Hydrogen  and  Air 


Probe 

Ft  From  ("enter 


1.833 

1.833 

1.833 

1.833 

1.833 

1.833 

1.833 

1.833 

1.833 

1.833 


Time 

Sec  After  Spark 

0.01442 
0.C1438 
0. 01490 
0.01476 
0.01488 
0.01476 
0.01451 
0.01780 
0.04146 


Probe 

Ft  From  Center 

1.666 

1.666 

1.666 

1.666 

1.666 

1.666 

1.666 

1.666 

1.666 

1.  666 


Time 

Sec  After  Spark 

0.01154 

0.01196 

0.01226 

0.01195 

0.01195 

0.01195 

0.01207 

0.01549 

0.03549 


Probe 

Ft  From  Center 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 

1.583 


Time 

Sec  After  Spark 

0.01106 

0.01148 

0.01202 

0.01171 

0.01146 

0.01171 

0.01146 

0.01500 

0.03037 

0. 0860 
0.0676 

0.  0601 

0.0617 

0.0798 


1.833 

0.1202 

1.666 

0.0913 

1.583 

] . 833 

0. 0812 

1.666 

0.0724 

1.583 

1.666 

0. 1331 

1.583 

0. 1243 

1.833 

0. 0693 

1.666 

0. 0626 

1.583 

1.833 

0.0693 

1.583 

0.  0626 

1.  833 

0.0760 

1.666 

0.  0672 

1.583 

1.  833 

0.0915 

1.  666 

0.0838 

1.583 

1.833 

0. 0851 

1.583 

0. 0759 

1.833 

0.0919 

x.  583 

0.  0789 

Table  A-5.  PVC  Pipe  Test  Conditions 

Pressure,  Pipe  Size, 


Run 

f/a 

nsia 

Temperature 

in. 

Fuel 

1 P 

0.064 

5.0 

Ambient 

2 

Propane 

2 P 

0.064 

14.0 

Ambient 

2 

Propane 

3 P 

0.  064 

14.7 

Ambient 

2 

Propane 

4 P 

0.064 

14.0 

Ambient 

2 

Propane 

5 P 

0.064 

14.7 

Ambient 

2 

Propane 

6 P 

0.064 

14.7 

Ambient 

2 

Propane 

7 P 

0.064 

10.0 

Ambient 

2 

Propane 

8 P 

0.064 

10.0 

Ambient 

2 

Propane 

9 P 

0.064 

10.0 

Ambient 

2 

Propane 

10  P 

J.  064 

10.0 

Ambient 

4 

Propane 

11  P 

0.  064 

5.0 

Ambient 

4 

Propane 

12  P 

0.064 

14.7 

Ambient 

4 

Propane 

13  P 

0.064 

10.0 

Ambient 

4 

Propane 

14  P 

0.064 

10.0 

Ambient 

4 

Propane 

15  P 

0.064 

14.0 

Ambient 

4 

Propane 

17  P 

0.064 

10.0 

Ambient 

6 

Propane 

18  P 

0.064 

14.7 

Ambient 

6 

Propane 

19  P 

0.064 

14.0 

Ambient 

6 

Propane 

20  P 

0.064 

5.0 

Ambient 

6 

Propane 

21  P 

0.064 

14.0 

Ambient 

6 

Propane 

22  P 

0.029 

14.7 

Ambient 

6 

Hydrogen 

23  P 

0.029 

14.7 

Ambient 

6 

Hydrogen 

24  P 

0.  029 

14.7 

Ambient 

4 

Hydrogen 

25  P 

0.029 

14.7 

Ambient 

2 

Hydrogen 

26  P 

0.029 

14.7 

Ambient 

2 

Hydrogen 

27  P 

0.029 

5.0 

Ambient 

2 

Hydrogen 

48 


Run 

Probe  1 

In.  From  Center 

Time 

Sec  After  Spark 

Probe  2 

In.  From  Center 

Time 

Sec  After  Spark 

Probe  2 

In.  From  Center 

1 P 

0.0 

0. 03679 

0.70 

0. 05873 

0.35 

2 P 

0.0 

0.01805 

0.70 

0. 02455 

0.35 

3 P 

0.0 

0.01847 

0.70 

0. 02369 

0.35 

4 P 

0.0 

0.01819 

0.70 

0.02223 

0.35 

5 P 

0.0 

0.01386 

0.70 

0.01795 

0.35 

6 P 

0.0 

0.01418 

0.70 

0.01848 

0.  35 

7 P 

0.0 

0.02133 

0.70 

0. 03458 

0.35 

8 P 

0.0 

0.02145 

0.70 

0. 03337 

0.35 

9 P 

0.0 

0.02110 

0.70 

0. 03512 

0.35 

10  P 

0.0 

0.03195 

1.60 

0.  04500 

0.80 

11  P 

0.0 

0. 02976 

1.60 

0. 08366 

0.80 

12  P 

0.0 

0.02451 

1.60 

0. 03146 

0.  80 

13  P 

0.0 

0.02953 

1.60 

0.  04419 

0.  80 

14  P 

0.0 

0.03085 

1.60 

0. 04329 

0.80 

15  P 

0.0 

0.02585 

1.60 

0. 03220 

0.80 

17  P 

0.0 

0.03866 

2.66 

0. 04500 

1.33 

18  P 

0.0 

0.03110 

2.66 

0. 03805 

1.33 

19  P 

0.0 

0. 02695 

2.66 

0.  04549 

1.33 

20  P 

0.0 

0. 05744 

2.66 

0.11134 

1.33 

21  P 

0.0 

0.02707 

2.66 

0. 06305 

1.33 

22  P 

0.0 

0.01024 

2.66 

0.01159 

1.33 

23  P 

0.0 

0.01200 

2.66 

0.  01679 

1.33 

24  P 

0.0 

0.00578 

i.60 

0.00741 

0.80 

25  P 

0.0 

0. 00378 

0.70 

0. 00518 

0.35 

26  P 

0.0 

0.00372 

0.70 

0.00512 

0.35 

27  P 

0.0 

0.00646 

0.70 

0. 01354 

0.35 

NOTE: 

Probe  1 - Located  1 ft  axially  from  spark  plug. 
Probe  2 - Located  2 ft  axially  from  spark  plug. 
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Table  A-6.  PVC  Pipe  Data 

L 


Probe  2 

In.  From  Center 

0. 35 
0.  35 
0.  35 
0.35 
0.35 
0.35 
0.35 
0.  35 
0.  35 
0.  80 
0.80 
0.80 
0.  80 
0.  80 
0.80 

1.33 

1.33 

1.33 

1.33 

1.33 

1.33 

1.33 
0.80 
0.  35 
0.35 
0.  35 


Time 

Sec  After  Spark 

0. 05867 
0.02444 
0. 02358 
0. 02211 
0. 01780 
0.01830 
0. 03434 
0.03277 
0.03488 
0. 04500 
0. 08427 
0.03061 
0.04442 
0.  04280 
0.03146 
0. 04305 
0. 03671 
0.04195 
0. 11085 
0. 05878 
0.01146 
0.01661 
0. 00729 
0. 00518 
0. 00506 
0.01354 


Probe  2 

In.  From  Center 

0.0 

C.O 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

0.0 

0.0 


Time 

See  After  Spark 

0.05861 
0. 02438 
0. ^2347 
0. 02211 
0.01773 
0.01824 
0. 03422 
0.03265 
0. 03463 

0.08439 
0.03049 
0. 04442 
0. 04268 
0.03134 


0. 11134 
0. 05829 


0. 00717 
0. 00518 


Probe  2 

In.  From  Center 

-0. 35 
-0. 35 
-0.  35 
-0.35 
-0.35 
-0.35 
-0. 35 
-0. 35 
-0.35 
-0.  80 
-0.80 
-0.80 
-0.  80 
-0.80 
-0.80 
-1.33 
-1.33 
-1.33 
-1.33 
-1.33 
-1.33 
-1.33 
-0.80 
-0.  35 
-0. 35 
-0.35 


Time 

After  Spark 

0. 05867 
0. 02438 
0. 02330 
0.02217 
0.01773 
0.01830 
0.  03422 
0.03277 
0.03488 
0. 04463 
0.08378 
0. 03061 
0. 04442 
0.  04293 
0.  03146 
0. 04305 
0. 03659 
0.  04171 
0.11134 
0. 05732 
0.01146 
0.01673 
0.00729 
0.00518 
0.00512 
0.01354 


I 


Time 

Sec  After  Spark 


Probe  2 

In.  From  Center 


Time 

Sec  After  Spark 


Probe  2 

In.  From  Center 


Time 

Sec  After  Spark 


APPENDIX  B 

CENTRIFUGE  COMBUSTION  PROCESS 


The  combustion  process  in  a dosed,  rotating  pipe  is  very  complex  to 
model.  The  model  must  include  flame  bubble  mechanics,  pressure  rise  due  to 
burning,  hot  gas  expansion,  and  constantly  changing  parameters  of  the  burned 
and  unburned  gases.  The  present  analysis  was  limited  to  the  initial  portion  of 
burning  and,  therefore,  certain  simplifying  assumptions  were  made.  These 
were: 

1.  The  temperature  of  the  burned  gas  was  constant 

2.  The  velocity  of  the  burned  gas  was  zero. 

After  ignition,  the  flamefront  travels  down  the  tube  at  an  observed  velocity, 
SQ.  The  cold,  unburned  gas  travels  at  S^.  caused  by  the  thermal  expansion  of 
the  hot  gas.  If  one  is  riding  the  flamefront,  the  cold  gas  is  entering  at  the 
bubble  velocity,  S^,  and  leaving  at  the  observed  velocity,  So>  From  the  con- 
servation of  mass,  neglecting  unsteady  terms,  we  get: 

P Sn  P. sS  Eq  (22) 

a B B o H ' 

where 

P cold  gas  density 
n 

true  flame  velocity 
Pj^  ~ hot,  burned  gas  density 

S = observed  flamespeed. 
o 

Neglecting  unsteady  terms  proved  to  be  a valid  assumption  in  this  case  because 
(iensity  changes  were  relatively  smaii. 

The  basic  bubble  equation  can  lie  written  as 


1 - 


P 


B 


9 

Ru)“  • 


Eq  (23) 


f>0 


where  nil  terms  hnve  been  defined  previously.  If  a value  for  A^C^/2  VR  is 
assumed,  the  differential  equation  can  be  solved  numerically  if  P^  and  p^  are 
calculated  as  functions  of  time.  Initially,  we  can  slater 


where 


p /R  T 
1 B/WB  B 

p/Rt" 

a a a 


P = density 
P = pressure 
T = temperature 
R = gas  constant 


Eq  (24) 


subscript  B bubble 
subscript  a cold. 


The  change  in  density  can  be  calculated  from  the  mass  fraction  burned 
during  each  time  iteration  assuming  constant  pressure  and  burned  gas  temper- 
ature. A new  pressure  is  calculated  from  the  fraction  of  mass  burned,  and  this 
process  is  repeated  down  the  tube.  A flow  diagram  is  presented  in  figure  B-l, 
followed  by  a listing  of  the  computer  program  and  a sample  result. 
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APPENDIX  C 

SWIRL  AUGMENTOR  MODEL 


The  swirl  augmentor  model  assumed  a hot  flame  bubble  was  driven  toward 
the  renter  by  centrifugal  forces  produced  by  the  swirling  flow  and  moved  axially 
by  the  flow  of  air  through  the  augmentor.  The  flame  bubble  location  was  calcu- 
lated from  the  axial  and  radial  movements  of  the  bubble.  The  two  movements 
were  assumed  to  be  independent.  The  axial  movement  of  the  flame  bubble  was 
calculated  from  the  axial  velocity  of  the  cold  unburned  gases,  which  increased 
down  the  duct  due  to  pressure  losses.  The  radial  movement  of  the  flame  bubble 
was  more  involved  and  will  be  discussed  in  detail.  The  model  allowed  a small 
fraction  of  combustible  gases  to  burn,  followed  by  an  expansion  of  the  burned 
gases.  Another  small  frar*ion  was  burned,  followed  by  another  expansion. 

This  process  was  repeated  until  all  of  the  combustible  mixture  was  burned. 

For  the  burning  step,  the  burning  rate  was  calculated  from  equation  4. 
This  required  three  terms  to  be  calculated  at  each  axial  station:  the  rotational 
rate  (u,),  the  density  ratio  (p^  ''p a ) . and  2VB'AE  Cp.  Of  these  three  terms,  the 
rotational  rate  required  the  most  analysis  and  will  be  discussed  first. 

To  calculate  o.,  an  analysis  of  the  swirling  flow  was  required.  The  flow 
immediately  downstream  of  the  vanes  can  be  analyzed  if  the  total  enthalpy  and 
the  entropy  are  assumed  radially  constant  and  the  radial  equilibrium  is  approxi- 
mated bv 

2 

f ^ = — Eq  (25) 

p dr  r / 

where 

g - Gravitational  acceleration 

/>  = Density 

p - Pressure 

V = Tangential  velocity 

r = Radial  position. 


From  Glassman  (Reference  6),  these  assumptions  result  in  the  following  equation 
downstream  of  the  vanes 


where  - axial  velocity. 


Eq  (26) 
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For  the  case  of  swirl  vanes  with  radially  constant  turning  angle  a , equation  26 
can  be  integrated  to 


Eq  (27) 


where  rc  represents  the  radius  of  half  the  total  cross-sectional  area  and  VUc  and 
VX(;  represent  the  tangential  and  axial  velocities  at  r^  Figure  C-l  shows  Vu/Vuc 
and  Vx/VXc  vs  radius  for  swirl  vanes  of  35  degrees.  For  the  case  of  small  vane 
angles  (<  35  degrees),  the  flow  can  be  approximated  by  radially  constant  tangential 
and  axial  velocities.  Thus,  the  angular  rotation  at  the  vane  exit,  w,  can  be  calcu- 
lated from  the  following  equation: 

V V Tan  a 

* = ~~  = “2— . Eq  (28) 


Figure  C-l.  Radial  Variation  of  Tangential  and  FD  89710 

Axial  Velocities  in  Scaled  Swirl 
Augmentor  Tested  Under  NASA 
Contract 


i 
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Downstream  ol  the  vanes,  the  cold,  unburned  gases  are  forced  toward  the 
center  due  to  the  expansion  of  the  burned  gases.  Since  angular  momentum  is  con 
served  along  a streamline,  V r - constant  along  a streamline.  The  rotational 
rate  downstream  of  the  vanes  was  calculated  using  this  fact,  and  the  analysis  is 
shown  below. 


m 


Jm 


Streamline 


Pilot  OflD? 


Burned  Gas 


Wl 

P2  A2  V2 


1 1 2 

ft;  (nRi 


2 2 
ff;  <ffR2 


- .7  RIV  ) M 


RIV")  M 


where: 


P Static  pressure 

R Gas  constant 

T Static  temperature 

M Mach  number 

i Ratio  of  specific  heats 

m Mass  How 

Rj,  R9,  and  RIV  from  diagram  above. 


II  T.  = '1'  :R.2-  -RIV2  A,  and  7 R 2 


- r RIV  = J 


Fro/ii  the  conservation  of  angular  momentum  along  a streamline: 


Substituting  equation  30  into  equation  31  we  get 


Eq  (31) 


A2  P2  MS 
P1  M1 


n RIV 


A2  + jtRIV 


Therefore, 


u)  = 
2 


V 


R, 


Eq  (32) 


Eq  (33) 


The  second  term  required  in  equation  4,  ^/a»  , was  assumed  constant 

throughout  the  burning.  If  the  fuel-air  mixture  was  ruch  that  the  mixture  burned 
at  a stoichiometric  flame  temperature,  the  density  ratio  was  0.393  for  an  inlet 
temperature  of  1044 °R. 

The  last  term  required  in  equation  4,  2VB/AB  CD,  was  calculated  from 
equation  13  assuming  a bubble  diameter  equal  to  the  pilot  dimension.  The  bubble 
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Reynolds  number  was  checked  to  assure  the  limit  size  was  not  reached.  The 
Reynolds  number  of  the  bubble,  assuming  a diameter  of  1. 15  inches,  was 


Rebub  " 174  SB 


Eq  (34) 


Therefore,  until  the  bubble  reached  a velocity  of  404  ft/sec,  the  limiting  size 

would  not  be  reached.  The  analytical  program  predicted  this  velocity  would  not 

be  reached.  Therefore,  2V-/A.,  was  assumed  constant  and  equal  to  0o  0296. 

ij  13  L) 

Once  the  terms  u>,  p„/p  , and  2VT3/AD  Cn  were  calculated  at  each  axial 

13  a Jo  13  u 

station,  the  radial  flamespreading  rate,  including  the  expansion  rate  term,  was 
calculated  step  by  step  in  the  following  manner.  Starting  at  time  zero,  the 
burning  rate  was  calculated  using  equation  4.  For  a small  increment  of  time, 
the  amount  of  flow  burned  was  calculated  and  the  weight  flowrate  of  the  hot  and 
cold  streams,  wh  and  wc,  were  adjusted  accordingly.  The  hot  gases  were  then 
allowed  to  expand  into  the  cold  stream.  Mathematically,  the  area  occupied  by  the 
cold  gas,  A , the  area  occupied  by  the  hot  gas,  A.,  and  the  static  pressure  were 
calculated  from  the  following  three  equations: 


i 

I 

I 
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where 


A * area 

w * weight  flowrate 

T^  = total  temperature 

Pt  = total  pressure 

P = static  pressure 

y = ratio  of  specific  heats 

R = gas  constant 

G = gravitational  acceleration 

Subscript  c = cold 
h = hot 
t = total . 


In  these  equations,  Pt  was  assumed  constant,  while  P^  deceased  axially  according 
to  the  fraction  of  mass  burned.  From  the  static  pressure,  the  axial  velocity  was 
calculated  and  used  to  determine  the  axial  movement  of  the  flame  bubble.  The  first 
step  of  the  burning-expansion  process  was  then  complete.  Each  succeeding  step 
was  done  in  the  same  manner  until  all  of  the  gases  were  burned. 

The  analytical  model  of  the  swirl  augmentor  was  incorporated  into  a com- 
puter program.  A listing  of  the  program  plus  theoretical  results  of  the  swirl 
augmentor  tested  under  NASA  Contract  NAS3-17348  are  presented  on  the  following 
pages. 
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00001170 

0017 

7f  * AT-AC-AH 

00001100 

0016 

KOUNT-l 

00001190 

0019 

P5TAT  * PSTAT1 

00001200 

0020 

KPT-0 

00001210 

0021 

1 K—I.O 

00001220 

0022 

X-0.0 

00001230 

0023 

IMPSTAT.CT.PTH)  PSTAT-PTH 

00001240 

0024 

PF«PS7AT/PTH 

00001250 

0025 

FH-0.0 

00001260 

0026 

CALL  IfNFLOCCAH»»A,HH.X.PAfX*XtFe*X.X) 

00001270 

0027 

AH  ■ 0.0 

00001260 

0028 

IF  IUH.GT .0.0.  ANP.FM.CT  .0.0) 

1 AH*VH»S0RT(«3.25*TH/37.17)/(PTH»FH)/I44.0 

00001290 

00001300 

0029 

MO- 1.0 

00001310 

0030 

X-0.0 

00001320 

0031 

pp-pstat/ptc 

0000133G 

0032 

FC-0.0 

00001340 

0033 

CALL  1 SNFLOIGAMMA ,NC«X*P?  ,X«X,Ff,X,X) 

00001350 

0034 

AC  * 0.0 

00001360 

0035 

IF(WC.GT.O.O)  AC«WC*<0BTI5^.n4:C/3?.17)/CPTC*FC)/l44.0 

00001370 

0036 

IFnnPT.fO.7)  PFTUPN 

00001360 

0037 

2-AT-AC-AH 

00001390 

0030 

CAl L SOLVIT (PSTATC, 70, P<TAT ,7 ,KfVNTf  If PT AVtO 1 ,SAVf02 1 

00001400 

0039 

KOHNT  « KPUfft  ♦ l 

00001410 

0040 

IFIKrUNT.CT.2000)  GO  TP  3 

00001420 

0041 

If  IPSTAT.lF.PTH.ANT.PSTAT.Gf  .PST/TD  CO  TO  1 

00001430 

0042 

4 WR1TM6,?)  WCtlf  .PTC.Af  tHr,WHfTHfPTM,*H*HH#AT,PSTAT.7 

00001440 

0043 

2 FORMAT  f • 1PRTF  IN  FUNCTION  P5TAT  - APfUHFNT*  APf  - •/ 

• • SUFI?.*/*  •••  P5TAT  ■ • * F 1? .7*  • AFFA  'PROP  * 

00001450 
* .fc  14.710000U60 

0044 

PF1UP0 

0000147C 

0045 

3 vr  tn<f  .4) 

0000)460 

0046 

4 '-c rmat ( • •••  nr  <rtUTioN  jn  pstatm 

0000)490 

0047 

cr  TO  5 

00001500 

f 

muc  UNITU)  ,TF  Afl  *INIT 

000G1MC 

C 

AT  1 

0C001520 

C 

tfaca  rv 

00001530 

0040 

too 

00001*40 

*1 

TT1 

MIN 

VANG 

GAMMA 

0 

TAOIA 

»1V 

AOV 

MU 

TT1 


EEf 


• 7. *00  IN 

« 164 A. 0 Off-K 

•>  0.11500 

■ 0.61000  A AT 

• 1.30000 

• 15.000  M 

« 4172.0  CIG-A 

• 2.000  IN 

« 6.333  IN 

• 30.000  MSI 

■ 16*4.0  OfC-P 


tm 

(SEC) 

0.0 

0.000500 

0.001000 

0.001500 

0.002000 

0.007500 

0.003000 

0.003500 

0.004000 

0.004700 

1.0 


IFNGTh 
If  T 1 
0.0 

0.1126 

0.2329 

0.36  t? 
0.M51 

0. 679? 
0 .6609 

1. C601 
1.2761 
1 .-666 


M AOlUS 
IE  T l 
C.527M 
0.49E* 
0.45*7 
0.41 0? 
0.3637 
0.3173 
0.2710 
0.2?*0 
0.1P01 
0 . 1667 


FtAMF  SPtft 
CFT/«FC  1 
0.0 

-47.E6M0 
-52.4669 
-57.4304 
-63. *271 
-70.200P 
-77.2003 
“M?.*  6 1 2 
“P5.0M35 
-M3. 776* 


rpir 


<FPfAK  V 
IM/StCI 
370.065? 
731.6797 
752.3521 
276.9634 
?0o#qP?? 
34*. 6943 
379.1725 
414.03*9 
445.94M2 
454.74Q0 


pure* 

I P Af'/SE  C 1 
490.0964 
530.4905 
601. 6f 31 
605.74*1 
614,7407 
062.1321 
1141.7440 
1349. 95M3 
1*33# If 74 
1551. fou 


JUl*  1IHF 

irtr.-Ri 

1644.0000 
1632.6671 
2090.0913 

2366.2(04 

?fcoo.64f? 

3037.0601 

3391.1707 

3747.2451 

4079.2623 

4177.0000 
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